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Theory of Diffraction in Microwave Interferometry 


JOURNAL OF RESEARCH of the National Bureau of Standards—B. Mathematics and Mathematical Physics 


Vol. 64B, No. 1, January-March 1960 


D. M. Kerns and E. S. Dayhoff ! 
(August 25, 1959) 


Microwave Michelson and Fabry-Perot interferometers are respectively considered as 
instances of: (1) A “reflection system’’, consisting of a radiating-receiving system and a 
reflecting object (e.g., a finite mirror); and (2) a “transmission system’’, consisting of a 
radiating system and a receiving system with an object (e.g., a Fabry-Perot etalon) inter- 
posed. The basic theoretical objective is the calculation of the amplitude and phase of the 
(time-harmonic) received signal in the systems considered. The electromagnetic field in 
space transmission paths is represented in terms of continuous angular spectra of vectorial 
plane waves, and the elements of the systems are described by means of suitable tensor 
scattering matrices (having both discrete and continuous indices). Needed scattering 
matrices are considered known; relationships to experimentally determinable data are 
outlined. The general case of either the reflection or transmission system is soluble formally 
in terms of a series of integrals stemming from the Liouville-Neumann series solution of 
certain integral equations. Formulas are obtained for models of the Michelson and Fabry- 
Perot instruments with arbitrary radiating and receiving characteristics. The theory and 
various features of the instruments considered, including Fresnel-region (or quasi-optical) 
behavior, are illustrated by means of examples obtained by choosing relatively simple and 


rather hypothetical analytical expressions for the radiating and receiving characteristics. 


1. Introduction 


Microwave versions of the Michelson and_ the 
Fabry-Perot interferometers are being studied and 
developed at the National Bureau of Standards with 
a view to high-accuracy measurements of wave- 
lengths, lengths, and the speed of light [1]. The 
data provided by these instruments are, however, 
subject to corrections for effects which are completely 
negligible in optical interferometry but significant 
in the microwave region. In this work—in contrast 
to the corresponding optical case—the principal 
such effect results from the wavelength not being 
negligibly small relative to the dimensions of the 
apparatus, and consequently the effects of diffraction 
upon the wavelength as observed by the instrument 
must be carefully considered in order to exploit 
fully the potential accuracy of the method. In 
this paper we consider these effects and present 
analytical tools suitable for making a correction to 
the apparert wavelength [2]. This should make 
possible the use of Fresnel region microwave inter- 
feometry for work of the highest accuracy. 

For the purposes of our analysis the pertinent 
aspects of the class of microwave Michelson inter- 
ferometers are represented by a “general reflection 
system.”’ Similarly, the microwave Fabry-Perot 
instruments are represented by a “general trans- 
mission system.’ (These basic arrangements are 
described more fully below.) Inasmuch as we wish 
to provide a theory inherently capable of dealing 
with high-accuracy experiments, a considerable 
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degree of generality is required to'avoid over-idealiza- 
tion. It is interesting that the required generality 
is such that the general reflection system could 
equally well represent a system consisting of radar 
and target, for example; similarly, the general 
transmission system po represent a point-to-point 
communication system. In view of the generality 
of the basic arrangements it seems likely that appli- 
cations of the theory will also be found in eas a 
other than the ones that motivated this work. 

The specific theoretical objective is the calculation 
of the amplitude and phase of the received signal in 
the systems considered. Theoretical expressions 
for the received signal will contain, at least implicitly, 
all available information regarding diffraction errors, 
intensities, etc. The analytical technique employed 
involves a transducer point-of-view, scattering- 
matrix formalism and constitutes a generalization 
of a part of the theory of waveguides and waveguide 
junctions. 

The arrangements to be considered involve 
radiating, receiving, and radiating-receiving systems 
as elements, which we shall occasionally refer to as 
“terminals” or as “‘terminal apparatus.’’ Such an 
apparatus is shown sliomndibellly in figure 1. The 
left-hand block in the figure may contain sources of 
power, detectors, and any other required auxiliary 
equipment, all of which is assumed to be shielded. 
This equipment is connected by means of waveguide 
to a waveguide-space transducer, represented by the 
triangular symbol in the figure. For simplicity in 
the formulation it will be assumed that the wave- 
guide is of a self-enclosed or shielded type. The 
arrangement and details of the apparatus are 
highly arbitrary, an essential but not very restrictive 
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FicurE 1. Schematic radiating-receiving system. 


condition being that the transducer proper possesses | 


the property of linearity with respect to electro- 


magnetic fields. (Details of the statement of 
a, pep are given in sec, 2, A.) 
n one of the basic arrangements considered—the 


general reflection system, illustrated in figure 2- 
radiation from a radiating-receiving system is re- 
flected or scattered by a general reflecting object at 
a variable distance d and is partly received (and 

artly reflected or scattered) by the same system. 

he Michelson interferometer can be considered in 
general terms as a reflection system in which the re- 
ceived signal is observed, as a function of d, by 
causing it to interfere with a suitable reference 
signal and detecting the resultant amplitude. 

In the laboratory version of the Michelson inter- 
ferometer being considered in connection with the 
preseat work, the radiation forms a_ well-defined 
beam (within its Fresnel zone), the reflecting object 
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Figure 2. Reflection system. 


If the scattering object is a finite mirror, for example, the arrangement represents 
amicrowave Michelson interferometer. 


is a large but finite mirror placed well within the 


Fresnel zone of the beam, and multiple reflections 
between the mirror and the radiating-receiving Svs. 
tem are relatively unimportant. It may be said 
that the instrument is operated in the Fresnel] op 
quasi-optical region. This mode of operation jm. 
plies the inequalities a>) and 
a is a measure of the cross section of the beam and 
\ is the free-space wavelength corresponding to the 
frequency of operation. Under these conditions the 
radiation in the space between radiator and mirror 
approximates a homogeneous, plane standing-waye 
and the received signal is approximately propor- 
tional to exp (2ikd), where k=2r/X. 


d<a*/y, where. 


The second basic arrangement considered is illys-" 


trated in figure 3. Here a structure or object that 
in general may reflect, transmit, and scatter is inter- 
posed between a radiating and a receiving system, 
If the interposed object is some form of Fabry-Perot 
“etalon” [3], the arrangement represents a general 
form of the Fabry-Perot interferometer. In this 


type of instrument the received signal exhibits am-* 


plitude variations sharply dependent on the spacing 
of the plates of the etalon, and this is attributable to 
interference among components of the radiation 
having experienced 0, 1, 2, reflections between 
the plates. 
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FiGcurRE 3. Transmission system 


If the scattering object is a suitable form of Fabry-Perot etalon, the arrangement 
represents @ microwave version of the Fabry-Perot interferometer. 


In the laboratory version of the Fabry-Perot be- 
ing considered in connection with the present work, 
the instrument is a Fresnel-region instrument in the 
sense explained above, and multiple reflections, ex- 
cept between the plates of the etalon, are again 
reiatively unimportant. 

In the theoretical approach used in the present 
paper, the (vectorial) electromagnetic field in a space 
transmission path is represented as the superposition 
of an angular distribution or spectrum of (vectorial) 
plane waves. Such a spectrum is in general deter- 
mined not only by the free-space radiation charac- 
teristics of a radiator, but also by the effects of all 
scattering elements involved (the radiator itself 
would be such an element, for example, in the case 
of a Michelson instrument with appreciable multiple 
reflections). Although the use of the plane-wave 
representation of the fields is suggested by the close 
approximation to a single plane wave that might 
exist in a Fresnel-region instrument, it is hardly 
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cessary to remark that in a practical case even the 
unmodified free-space radiation spectrum is likely to 
be extremely complicated in detail. This is because 
of such things as edge effects, geometrical imperfec- 
tions, and limitations in the design and construction 
of lenses and horns. 

The concept of the plane-wave resolution of the 
field permits one to form a useful qualitative picture 
of the origin of an effective wavelength. Consider, 
for example, a Michelson instrument. Each ele- 
mentary wave in the spectrum whose normal makes 
an angle @ with the line along which translation of 
the mirror is measured has an effective wavelength 
rh sec 6; the resultant effective wavelength can be 
thought of as a kind of average of the elementary 
contributions over the existing spectrum. Thus one 
might anticipate, for example, that the resultant 
effective wavelength should rather generally tend to 
be greater than the free-space wavelength. 

A main section of this paper, section 2, is devoted 
to establishing a seattering matrix formalism for 
the description of radiating-receiving systems of the 
general type considered. In this connection it 
helpful —at least to one acquainted with the theory 
of waveguide junctions—-to regard the space side 
of the waveguide-space transducer as a waveguide 
of infinite cross section. A scattering matrix having 
both diserete and continuous indices, which corre- 
spond respectively to the mode in the ordinary 
waveguide and the continuum of modes in the wave- 
guide of infinite cross section, is required. The 
expression of reciprocity, which seems to be of in- 
terest in itself in the class of problems considered, 
is established. 

For the purpose of the present work, the scattering 
matrices of the radiating-receiving systems involved 
are considered known; the determination of the nec- 
essary data for a practical application is considered 
to be an independent problem, experimental or pos- 
sibly theoretical. (Fortunately in the practical cases 
under consideration, seattering by the radiating- 
receiving systems involved is of minor importance, 
and it thus appears that the necessary data can be 
obtained from suitable measurements directly or 
indirectly determining the far-field radiation charac- 
teristics. Such measurements would vield the neces- 
sary information about the underlying diffraction 
problem, however complicated it might be.) 

In section 3 the seattering matrix formalism is 
extended and applied to the calculation of the re- 
ceived signal in the basic arrangements described 
above. The general problem involving multiple 
reflections is soluble formally (in the sense that the 
problem is reducible to quadratures) in the form of 
a series of integrals, which stem from the Liouville- 
Neumann series solution of a system of integral 
equations. Formulas are derived for models of the 
microwave Michelson and Fabry-Perot interferom- 
eters in which the radiating and receiving spectra 
are arbitrary. Finally, these formulas are illustrated 
by means of examples obtained by choosing specific 
and relatively simple analytical expressions for the 
spectra involved. 


ne 


is 


2. Scattering Matrix Description of Radiat- 
ing-Receiving Systems 


A. Basie hypotheses. To the general description 
of the radiating-receiving systems given in section 1 
we add the following details. 

We choose a terminal surface Sp in the waveguide 
and a supplementary surface Sj (such that S)+S; 
forms a closed surface) coinciding with the shielding 
around source, receiver, ete. (See fig. 1.) We 
further choose a rectangular coordinate system 
Oxyz such that coordinate surface z=0, which we 
denote by S,, will serve as a terminal surface on the 
space side of the waveguide-space transducer. 
As a convenient artifice we employ an_ infinite 
hemisphere S, lying in 2<0 and centered at O. 
The interior—.e., the domain of the electromagnetic 
field—of the transducer is the region V bounded 
externally by S,+S, and internally by Sy+Sj. In 
figure 1 the structure of the transducer considered 
lies in V; other structures or objects may also be in 
\ but no attempt is made to illustrate this possibility. 

The whole of the space and structure within V is 
counted simply as a (decidedly) nonhomogeneous 
medium which may also be dissipative and aniso- 
tropic. Linearity assumed and is essential; 
reciprocity, in the sense that tensors describing the 
medium within V are required to be symmetric, is 
assumed and is useful but not essential. 

In the space transmission paths (i.e., in 2>0) the 
medium is to be homogeneous, isotropic, and non- 
dissipative, as well as linear. For the time being it 
is assumed that these properties hold for arbitrarily 
large 2; the interposition of elements (such as a 
reflecting object) will be considered later. (The 
theory is not actually restricted to transmission 
media having the ideal characteristics listed, since, 
as will be seen, an arbitrary linear transmission 
medium is to be treated as a suitable interposed 
element. ) 

It is assumed that the electromagnetic field quan- 
tities vary harmonically with time ¢ at frequency 
w/(27). We employ the usual complex electric and 
magnetic field vectors, E, H, which are functions of 
the position vector r of Oxyz, and omit the time 
dependent factor exp (—twt). 

It should be observed that in taking S, to be 
hemispherical it is tacitly assumed that the problem 
is not two-dimensional (by a “two-dimensional” 
problem is meant one in which all quantities are 
independent of one rectangular coordinate perpen- 
dicular to the z-direction). To avoid unduly com- 
plicating the discussion, two-dimensional problems 
are not considered in the main part of this paper. 
However, two of the examples at the end of the 
paper are two-dimensional. For convenience in 
discussing these and other two-dimensional cases, 
some of the more important or less obvious of the 
needed formulas are summarized in an appendix 
(sec. 4, B). 

B. Representation of field on So. It, as is assumed, 
only one mode (a propagated mode) is of importance 


Is 








in the waveguide in the neighborhood of So, then 
from waveguide theory it is known that the trans- 
verse components E,, H, of E, H on Spy may be 
written 


E (r) = (ao +5) e(r), 
H (r),= (ay— bp) h (r), 


(ron S,) (1) 


where e, h are real basis-fields for the mode involved. 
The basis-fields are subject to the impedance normal- 
ization 


h (r) =n Xe (r) (2) 
and the power normalization 
J e(r) <h(r)-ndS—4r' no, (3) 
Sy 
where n is the unit normal on S, drawn into Vo and 


nm is the wave-admittance for the mode involved. 
These equations implicitly define the quantities 
a>, bo; it can be verified that ad», by so defined are 
respectively linear measures of the electric field of the 
incident and emergent traveling-wave components 
of the waveguide field at Sy). The net time-average 
power input 7) at S) is given by 


Py=Re| LEX H-ndS=2r? no (\ao|?—|bol2), (4) 
So 


where Re denotes that the real part is to be taken 
and the superposed bar denotes the complex con- 
jugate. 

C. Representation of the field in the region z>0. 
As mentioned in the introduction, the electromag- 
netic field in the region z>0 is to be represented as a 
superposition of plane-wave solutions of Maxwell’s 
equations. This type of representation well 
known, at least in the case of solutions of the scalar 
wave equation [4]; the generalization to the vector 
electromagnetic field offers no particular formal 
difficulty. 

The electromagnetic field in the region under 
consideration satisfies Maxwell’s equations in the 
form 


Is 


Vv < E=ieuH, v< H -iweE, (5) 
where uv, € are constant real scalars representing 
respectively the permeability and the permittivity 
of the medium. (Rationalized MKS units are 
employed.) We derive our basis fields from the 
general plane wave 
E=T exp (ik-r), 
(6) 
H= (wu) ~'k XT exp (7k-r), 


which is a solution of (5) for any propagation vector 


| 


k such that k?=w*ye and any vector T satisfying 





the transversality relation k-T=0. (In spite of 
this occurrence of “transversality,” in what follows 
the term ‘‘transverse’”’ will always mean transverse 
with respect to the 2-direction.) 

The propagation vector will be regarded as q 
function of its transverse components k,, ky; the 
z-component is then 


k,=+7, (7a) 

) 2 9\1 . ° 
where y=(k?—k?—k?)?. 1t will be convenient to 
denote the transverse part of the propagation 


vector by K, so that K=s,e,+-k,e, and 


y= (k?—K?)}, (7b) 
Since k,, k, must be allowed to vary independently 

in the range (—o, o), real and imaginary values 

of y will occur. y will be taken positive for K?<k. 

positive imaginary for A*>k®. Superscripts “+” 

and ‘*—”’ will be used when it is desired to indicate 

the choice of sign associated with k,. 

In virtue of the relation k-T=0, (6) vields just 
two linearly independent fields, hence just two basis 
fields, for any given k. The appropriate polariza- 
tions for the basis fields are those with the electric 
vectors parallel or perpendicular to the plane of 
k and e., which is the plane of incidence for a ray 
incident on any plane z=const. This choice of 
polarizations corresponds to the preetice in elee- 
tromagnetic theory in deriving the Fresnel equa- 
tions, for example; it also corresponds to the resolu- 
tion into “transverse magnetic” and “transverse 
electric”? modes of waveguide theory. 

In setting up the desired basis fields it is convenient 
to employ the transverse unit vectors 

x, —K/A, x. e.k,, (8) 
which are respectively parallel and perpendicular 
to the plane of kande,. (The notation is illustrated 
in fig. 4.) As a temporary abbreviation we put 
u=—exp (/k*-r). For the “&)” (or TM) compo- 
nents we take T=«,+ Ay 'e. and obtain from (6) 


E* =|«,-Ky~'e,]u*, 
(9) 
H;? me, XK\u~, 
where n,;=we/y. For the “/,” (or TE) components 


we take T=«, and obtain from (6) 


Es = «,u* 
(10) 
HF =[+ me, < «,+ K (wy) ~'e,}u* 
where m.=y/(wu). Among other similarities it may 


be observed that 7, n. are ‘‘wave-admittances” that 
correspond exactly to the wave-admittances en- 
countered in the theory of waveguides with discrete 
modes. The expressions in (9) and (10) are essen- 
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Figure 4. Illustrating k, K, x, and xe. 
tially the desired basis fields (2-components, in- 


cluded here, will be dropped later). The normaliza- 
tion and orthogonality properties of these fields are 


of course implic it in the expressions themselves. 


introduce spectral density functions 
a,,(K) for outgoing and incoming 
and form the general super- 


We now 
b,,=6,,(K) and a,, 
waves, respectively, 
position 
. > 
E(r) = >> (6, Ex +-4,Ex )dK, 


> (11) 





H(r)— | 33 (b His +-d,Ha)dK. 
4 


Here and subsequently in expressions of this type 
summation over the two values of the polarization 
index m and integration over the domain of the two 
variables k,, k, will be understood. The E, H given 
by (11) will satisfy Maxwell’s equations provided 
merely that the necessary differentiations can be 
taken under the integral sign. However, for our 
purposes, the z-components of the fields are redun- 
dant. By discarding the z-components we obtain 
the much more convenient and explicit expressions 
for the transverse field components 





« 
E(r), | >> (6,¢6'7 -a,,e~ *?*)«,,e'® RdK, 
> (12) 
H(r), | 2 ae ™ mn?) 1n@eX Kye ™ FdK, 
. J 
where R denotes the transverse part of r. These 


equations exhibit E,, H, as two-dimensional Fourier 


transforms of the quantities multiplied by exp 
((K-R) in the respective integrands; from the inverse 
transformations one may obtain 


iyz 
b,, (K) - is fi (E+,'HxXe,)e"® BIR, | 
, (13) 


an(K) = 3m: |(E- ng 'HXe.)e~™® BAR. 





J 





| 
| 


| the transducer as a whole. 








Here and subsequently integrations with the differ- 
ential symbol dR are to be taken over all values of 
zand y for some fixed z2>0. These equations enable 
one in principle to evaluate a,, 6, from a given 
distribution of E,, H, on any transverse plane. 

The case where there is no field incident on S, from 
the “right” [i.e., where a,,(K)=0] is of particular 
interest in what follows, and we shall list several 
relations holding in this case. From (12) and (13) 


E(r), [BaKemuK, (14) 


B(K) = (2r)~*e * {E().¢ iK-RdIR (15) 


where B=B(K) is the transverse vector }5,«;+ box. 
The time-average power, P,, radiated into the half- 
space z>0 will also be of interest. This power is 
the same as the time-average power flux in the 
+-2-direction across the surface S,. Hence, from 
(12) [with a,,(K) =0], 


JKR 


Re |E <H-e.dR=2r° ST (K) | by (IK) |20K, 

(16) 
indicated by the notation below the 
evanescent waves are excluded from 
the integration. This equation is of a type known 
as Parseval’s formulas; it may be derived formally 
with the aid of the rule that 


where, as 
integral sign, 


| e' *2—*)) 4d 7 =2nb(k,—k) (17) 


where 6 denotes the so-called impulse or 6-func- 
tion [5]. 

D. Definition of scattering matrix. Having set up 
representations for the fields on the two terminal 
surfaces of the waveguide-space transducer under 
consideration, we are now in a position to consider 
It may be assumed that 
a set of out-going wave-amplitudes [by and 5,,(K)] 
will be determined by a set of incident wave-am- 
plitudes [a> and d»(K)]. In fact, since the electro- 
magnetic system under consideration is by hypoth- 
esis a linear system, the relation between the set 
of out-going wave-amplitudes and the set of incident 
amplitudes must be a linear relation. We write this 


b= Sette | 33Sy,(m,K)a,(K)dK (18a) 


b,,, (IK) Syo(m K)ay+ [ SSi(m,Kyn,L)a, (Lal, 
(18b) 


thereby defining the scattering matrix for the trans- 
ducer considered. Figure 1 may be helpful i in fixing 
the significance of the quantities involved in (18). 
It is convenient and it seems appropriate to use the 








term ‘“‘matrix’”’ here even though one must think of 
rows and columns labeled both by discrete indices 
and by indices having continuous ranges. Evi- 
dently the functions Sy, Si, and Sj), respectively 
embody the receiving properties, the radiating 
properties, and the (space-side) scattering properties 
of the transducer involved; Soo is an ordinary wave- 
guide reflection coefficient defined at the terminal 
surface Sp and expressing the ‘antenna mismatch.” 

Equation (18) can be represented in terms of Soo 
and linear functional operators Sp;, Sj, and S),; such 
that 


bo = Sooo + Sod, (19a) 


A A 


b=Sypao+ Syd, (19b) 
where 6 and a@ are understood as function vectors 
corresponding to the functions 6,,(K) and a,,(K). 
This compact notation is used later primarily as a 
convenience in some of the more formal and general 
parts of the discussion. 

EK. Reciprocity. Radiating and receiving charac- 
teristics are related by the reciprocity condition, 
which here takes the form 

—noSo(m, K)=n,,(K)Sy(m, —K). (20) 
Since this particular form of the condition appears 
to be new, a derivation is given below (appendix). 
The occurrence of the factors No and nm(K) in (20) 
can be regarded as a consequence of the particular 
normalizations adopted in setting up the basis fields. 
(The same type of relationship holds for the elements 
of the scattering matrix of an ordinary waveguide 
junction with discrete modes [6].) The occurrence 
of the argument —K in one side of (20) means that 
the equation relates radiating and receiving charac- 
teristics in the line of a given propagation vector k 
(if a radiated wave has propagation vector k, 
the received wave in the same line has the propaga- 
tion vector —k). 

The plane-wave into plane-wave scattering func- 
tion Sy, is also subject to reciprocity ; the relation is 
nm(K)S),(m,—K; nL) =7,(L)S\,(n,—L: mK). (21) 
Comments similar to those following (20) apply here 
also. This relation is not used in the present paper 
but it may well be of interest in other diffraction 
problems. 

KF. Determination of scattering matrix. Although 
the basic viewpoint of this paper is that Soo, So, ete., 
are to be considered known, it seems well to take some 
note of the problem of determining these quantities 
from empirical data in the arrangements of particular 
Interest. 

Concerning So and S;, there is not much to be 
said. So is not only relatively easily measurable 
but also experimentally controllable (by means of 
tuning elements), whereas in general the same is not 
at all true of S,;,. However, the desired operating 
condition that S,, be effectively negligible is ap- 


| 


} 
proximately attainable and to some extent subjeg 
to experimental verification (by observation of thp 
effects of multiple reflections). 

Concerning So, and Sj we first note: (1) Eithe 
one of these functions may be determined readily 
from the other with the aid of the reciprocity relg. 
tion; (2) in the present context these functions need! 
be evaluated only for A?’—/*, the effects of evap. 
escent waves being avoided by keeping d >» (this 
is not inconsistent with Fresnel-region operation, 
cf see. 1). 

A direct approach to the determination of §) 
is implied by the definition (18): So, represents the 
received signal by as a function of the direction and” 
the polarization of incident plane waves of suitably 
normalized amplitudes. 

According to (18), Sio(mK) represents — the 
transverse components of the vector spectrum of 
outgoing waves under the conditions a,,(K)=0 and 


ady=1. If we define the transverse vectorial spectrum 
S,o(K) = Sw, K)«, + Sy,(2,K) x, (29 
then, from (15), 
3 (K | {pe K-R 
S,,(K) E,(Ri« : dR: (9°3) 
4rd, 


that is, Sie may be represented as the Fourier trans 
form of the transverse components of the electric 
field obtaining on the reference plane in the absence 
of incident waves, normalized to unit @. This is, 
of course, essentially a vector form of the well-known 
relation between “aperture” distribution and _ spee- 
trum. 

The vectorial spectrum Sj, (K) is also closely , 
related to the far electric field by a well-known 
type of relation. Under certain restrictions E(r), 
has for large r the asvmptotie form [7] 

(lf), ewes 2rik cos 6 B(RE re" /r; (24) 
the angle #@ introduced here is the polar angle ofr 
relative to the z-axis. By rewriting this equation 
and dividing by d) we obtain a formula for S,(K) 
in terms of the asymptotic form of E-: 

95) 


«-) 


S,)(K) =71(2ra,) ~'y~'r e~ *E(Kr/k) ¢ asymp > 
for A? <k?. 

Finally we note that the familiar “power radiation 
pattern” or “polar diagram’’ of antenna theory, 
defined as radiated power per unit solid angle as a 
function of direction, is given by 

1 Fa = an 

p=2(e/u)? (rk cos 6)? b-b. (26) 

where b=B-+-b.e. and 4, 

plete vectorial angular spectrum, including the 

z-component). Clearly this equation is not suf 

ficient by itself to determine S,)(K); polarization 
and phase information is required in addition. 


-B-K/y (b is the com- | 


fir 


eq 
ob 


Is 


tr 


(t 
ca 
or 


”" 


subjec 
of the 


Kithe 
readily 
Vv rela. 


IS need! 


evan. 
X (this 
Tation, 


of Sy 
its the 


on and * 


lit ably 


Ss the 
um of 
0 and 


ctrum 


(99 


trans- 
lectric 
SENCE 
ls 1s, 
hownh 
spec- 


losely 
nown 


Er), 
(24) 
ofr 


ation 


ro K) 


(25) 


ition 
Ory, 
as a 


3. Applications 


A. Reflection systems. One obtains a form of the 
first basic arrangement described in section 1 by 
placing an infinite plane reflecting surface “in front 
of” a radiating-receiving system of the type con- 
Ghaed in section 2. This represents a problem 
of intermediate complexity, from which the basic 
equations for the Michelson arrangement: may be 
obtained by specialization. If the reflecting surface 
is at zd and has reflection coefficient p(m, K), then, 
transforming the reflection coefficient to the plane 
-—0, we have 


a,(K) —p(m,K)e%®5, (K) 
(the dependence of y on its arguments is now indi- 
eated explicitly). Upon substituting (27) into (18) 
one obtains 


bo = Sooo 4 | > Sor (mK) p(m ,K) e?'"4b,, (KR) dK, 


(28a) 


b(K) = Syy(m, Kay 


+ {= Si. (m.K3n,L) p(n, Lye? (Lyd. (28b) 
The last line represents two simultaneous, inhomo- 
geneous, linear integral equations for the determina- 
tion of b,,(K) (do being prescribed). The Liouville- 
Neumann series solution of these equations may be 
obtained by a process of successive approximations. 
For the first approximation one takes 

by, (KR) = Syo(m,K) ap; (29) 
the second approximation is obtained by substituting 
the first into the right-hand side of (28b), 


b? (K) = Sio(m, Kj ay 


bdso(n,L)dL; 
(30) 


Fedo Si(m, K;n, L) p(n, L) e* 


and so on: 6% (K) accounts for the first n reflections 
at the reflecting surface. Once 6,,(K) is obtained, 
approximately or otherwise, it is to be substituted 
into (28a), thus determining the received wave- 
amplitude 6) in the waveguide at So. 

Useful approximate equations describing the be- 
havior of the Michelson instrument may now be 
obtained. The appropriate conditions are (1) that 
the effects of multiple reflections be negligible and 
(2) that the reflecting surface at z2=d be a mirror, 
for which we may put p=—1. The first condition 
means that (29) is already a good approximation 
for b,,(K); this substituted into (28a) yields 


b, r — , r 9 : . y - 
Soo ee Sy (m Kye>o™48,.(m, K)dK. (3 


dy e 


1) 


7 


The second term on the right will be called the “re- 
flection integral’ and denoted by #(d). Either So, 


or Sip may be eliminated by means of the reciprocity 
condition (20); for the purposes of the present dis- 
cussion it seems preferable to eliminate Sy. We 
then have 


® (d) {x Om (K) Sip (mK) Sy (m,—K) 2 MU K, 
(32) 
This is the main equation for the Michelson. It 


will be illustrated below by means of examples ob- 
tained by choosing specific mathematical expres- 
sions for Sio(m,K). 

Problems involving reflecting or scattering objects 
other than an effectively infinite reflecting surface 
are important not only in microwave interferometry 
but also in other fields. Consequently the following 
formulation of the general case where the plane 
reflecting surface considered above is replaced by an 
arbitrary reflecting object is of interest. Let the 
general reflecting object be characterized by means 
of a scattering function R(m,K; n, L) defined with 
>=0 as terminal surface. Instead of (27) we now 
have the linear transformation 


a, (K) = | S72 R(m,K;n,L)b,(L)dL (33) 
as boundary condition, In operator notation (33) 
is written 
a lo 
and the equations corresponding to (28) are 
bo =S9@)4 S,, Rb, (34a) 
h S odo T Sy 2b. (34b) 
The Liouville-Neumann series solution of (34b) 
may be written 
b= SyAo+ LSA, (35) 
where 
L=>> (8h) (36) 
p=1 


is the operator corresponding to what is called the 
“resolvent kernel” in the theory of integral equations, 
Finally, for b) we obtain 

So RS oy + So RLS Ay. (37) 


b. = Sudo 


In this form the last term (specifically the operator 


L) accounts for multiple reflections. 

An instance of (37) is the case of a large but 
finite mirror in the Michelson interferometer, which 
has been considered using an approximate expres- 
sion for PR and neglecting multiple reflections [8]. 








B. Transmission systems; Fabry-Perot interferom- 
eter. We now consider a general transmission sys- 
tem consisting of a radiating system and a receiving 
system with an arbitrary (electromagnetically linear) 
intervening structure or medium. For the active 


terminal the pertinent descriptive equation is 
(19b), which for convenience is repeated here: 
b= Soo +S,,4. 195) 


This equation is understood to be set up with 

reference to a coordinate system Oxyz and space- 

side reference plane z=0, as detailed in section 2. 

Using primes to distinguish quantities associated 

with the passive terminal, we may write 
b= Sua’, 


(39a) 


b’—S,,a’, (39b) 
as the equations corresponding to (19). For these 
equations the space-side reference plane is z=d in 
the above-mentioned coordinate system and _ the 
eneral arrangement is shown in figure 3. Next, 
fet the structure and/or medium between the 
terminals be described by a set of linear operators 


T;;, defined with respect to z=0 and z 
planes, such that 


=d as reference 


yy one 
T Tob ’ 





(40) 


| we have 


@’ =Trb64+ Tb’. 


(The Fabry-Perot interferometer considered below 
will furnish an example of these equations. It will 
be a very special example, however, since specular 
reflection and transmission will be assumed, so that 
the operators 7’, will be diagonal and (40) is then 
reducible to a family of ordinary equations.) 

A method of solving the problem described in 
(19), (89), and (40) may be indicated as follows. 


The result of eliminating 6 and 6’ from (19b), (39b), 
and (40) may be put in the form 


a 


Ti, Soo TS TS | @ 
= (41) 
a’ Tx, Soo Ty Si TeSi tla’ 


As this form suggests, these equations may be solved 
for @ and 4’ by a process of successive approxima- 
tions similar to that used above. The received 
wave amplitude bj is directly determinable from 
(39a) as soon as 4’ is known. (by is also directly 
determinable.) 

To pass to the consideration of a highly simplified 
model of the Fabry-Perot interferometer, we first 
assume that S,,; and Sj, are effectively so small 
that reflections’ at the terminals of the system may 
be neglected; that is, we assume that the second 
term on the right in (41) may be neglected. We 





ality, since the planes 2 


then have the explicit expression 


A a 


191 Soo, 


a’ (42) 


7 


and consequently for the received wave-amplitude 
we have 


bo= Soi Tr Si0o. (43) 
Next, we assume that the structure described by 
the 7;,; is a Fabry-Perot “etalon”, consisting of g 
pair of elements corresponding to the two plates of 
an optical Fabry-Perot interferometer. Such ele. 5 
ments might be, for example, perforated metal sheets 
or stacked quartz plates [8]. We assume that each 
element is symmetric with respect to the z-direction 
and characterizable by means of a (specular) reflec. 
tion coefficient p(m,K) and a (specular) transmission 
coefficient 7(m,K) defined at the symmetry plane 
of the element as reference plane. It is not as., 
sumed that the elements are lossless; p and rf are 
subject merely to realizability conditions for passive 
elements. We let the two elements be located go 
that their reference planes coincide with z=0 and 
with z=d, respectively. It may be noted that so 
locating the elements implies no real loss of gener- 
0 and z=d are arbitrarily 
located with respect to the physical arrangements 
with which they are respectively associated. From 
the symmetry of the etalon as a whole with respect 


to z=d/2 and the fact that the 7;; must be diagonal 


7 Y b,d(K L)?,,(m,K) (44a) 
T= T2= 8 nnd (K—L)t.,(m,K), (44b) 
where 6(K—L) means 6(4,—/,) 6(k,—/,) and 
1+ (7? oe dl 
ti=p | — peri ’ (45a) 
2,ivd 

to a (45b) 

1— p*e*'™ 


as may be found by ordinary methods. In these 
expressions y, p, and +r in general depend upon K. 
When evaluated for K=0 and simplified somewhat 
as they may be for lossless elements, the expressions 
reduce to ones frequently used in discussions of the 
optical Fabry-Perot interferometer. 

Using (44b) and writing out (43) we obtain 


Neue, | 52 8%, (m,K) t, (mK) Syo(m,K)dK; 


finally, inserting (45b) and defining ¥(d)=b)/ay, 


V(d)— (= Sou ( m,K)r°(m,K)e'™®4Sio(m,K) dK 


] —p?(m,K) era 
(46) 


V 


(42) 
plitude 


(43) 
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This is the ‘‘transmission integral” for a Fabry- 
Perot interferometer. An example of the analytical 
evaluation of (46) will be given below, assuming 
constant p and 7 and choosing very simple expressions 
for So, and Sy. (It may be observed that in regard- 
ing (46) as a function of d, it is implicit that the 
terminal apparatuses remain fixed relative to the 
reference planes with which they are respectively 
associated. ) 

The examples that follow have been chosen to 
illustrate various features of the theory and of the 
instruments considered. Inasmuch as the examples 
are rather hypothetical in nature and mainly illus- 
trative, no thorough or rigorous discussions are 
attempted. 

C. Examples, Michelson case. 

C.1. Dipole. The following example seems well- 
suited to illustrate the theory, inasmuch as it 
involves both TE and TM field-components in a 
fairly complicated way, the integrals involved can 
be evaluated, and the form of the answer can be 
anticipated. In this example the radiated field is 
assumed to be identical to that of an elementary 
electric dipole of moment p located at O. 

To find the angular spectrum of the electric field 
we may proceed as follows. The appropriate Hertz 
potential is [9] 


TI = (dre) ~'pe“*’/r; 


the representation of the spherical wave exp(ikr)/r 


in terms of plane waves is [10], for 2>0, 
e*" r= —(2ri)- ‘| _~ leik-rdK; 
where k=k* is understood. Since E=VXV <I 


we have 


E=0] kx (kxp)y te dK, (47) 
where ('=(Sx*ei)~'. Hence the complete vectorial 
spectrum is b=Ck&(kxXp)y"'. This result holds for 
all K; immediate confirmation for A?<k? may be 
obtained from (25) using the asymptotic form 
—k?(4are) ~' X (rX p)e™ /r* 

for E [9]. 

According to the definition of Sj(m, K), we have 
in this example ‘ 

Syo(m, K) = (C/do) &m + [KX (KX p) ly? (48) 

for m=1, 2. These spectral components are to be 
substituted into the reflection integral (32) for the 
Michelson instrument. The coefficient of exp (2i7d) 
in the integrand of (32), after some vector-algebraic 
labor, is found in the present instance to be expres- 
sible in the form 


—we(C'/ao)? p- [KX (KX p’) |y7?, 


532245—60— 2 


| where C; is independent of d. 





| 





, 


is the negative mirror image of p, with 


where p 
—Pz,—Py, P:). Thus for (82) we have 


components ( 
© (d) = —weny ‘(C/a)? p fx < (kX<p’)y~'e"dK. (49) 


Comparing this expression with (47) it is seen that 


Cwe 
nos, 





$(d)=— p- E’\0, 0, 2d), (50) 


where E’ is the electric field of a dipole of moment 
p’ located at O, and E’ is evaluated at the image of 
O in the reflecting surface of the instrument. This 
result is indeed of a form that might be anticipated. 

If for simplicity we take p to be transverse, (50) 
becomes more explicitly [9], 

&(d) =C;,| (2ikd) ~*— (2ikd) ~? +- (2ikd) —Je?™, 

It is clear that the 
result in this example is not physically meaningful 
as kd—»0, for the magnitude of ® can not properly 
exceed unity. This defect is attributable to the 
neglect of scattering or re-radiation by the dipole; a 
plausible extension of the theory of this example 
taking scattering into account gives results qualita- 
tively well-behaved for all values of kd. 

C.2 Two-dimensional Gaussian.—To provide a 
reasonably simple analytical illustration of the be- 
havior of a Michelson instrument operated in the 
Fresnel region, we consider a two-dimensional case 
with the pure TM “Gaussian” spectrum 


bs (k2) =A PH, 
(51) 


The subscripts 1, 2 here refer to e,, e,, respectively, 
and A is an arbitrary amplitude. The correspond- 
ing distribution of E, on the reference plane z= 
is also Gaussian, 


b.(k,) =(). 


’ 2; (902 
E,=Const. e~*/@, 
E,=0, 


as follows from (14,) (see appendix, B). ‘The param- 
eter a is, in a well-known manner, a measure of the 
sharpness of the angular spectrum and a measure 
of the width of the distribution of E,. 

In this example it will be interesting to determine 
the absolute magnitude of © explicitly. We need 
the relation between |ay| and |A| and we obtain this 
from a consideration of energy balance under free- 
space radiation conditions. If the fraction A of the 
net input power at S, is radiated, we have 


k 
TN \A\"(1 = |Soo|”)h =r{ > 2 1m|Om(Kz) \*dkz, 
-k 
(52) 








where the left- and the right-hand sides of this equa- 
tion come from (4,) and (16,), respectively. 
The appropriate two-dimensional form of (32) is 


(dd) =—No {Salk So m, kz) Sio( m,—k,z)e aid “dk, 


(53) 


. 


where we still have n,;=—wey™ 
the present case this becomes 


and No= Y( wp) 


“p-—2™ 
Nolo 


fr’ exp (—a*k? + 2iyd)dk,. 


Letting ¢ denote the phase of A/a) and using (52) to 
eliminate |A/a | from the last equation, one obtains 


fr lexp (—a°’k?-4 
| Soo|?)h — 


*k 
| y- exp (—a’k2)dk, 
e k 


Qiyd) dle, 
&(d) = 


e*'4(] 


(54) 


If ka>>1, this expression yields (0)~—e*'¢ 
(1—|Soo\A—a quantity that may approximate unity 
in magnitude. 

Since we are interested in Fresnel-region behavior, 
an asymptotic expansion of the numerator of (54) 
in terms of inverse powers of a is appropriate. This 
asymptotic expansion may be found with the aid of 
Watson’s lemma, as given by Jeffreys and Jeffreys 
[11]. One finds 


1—2ikd 
©(d) = C,e2**4] 1 4+-— = 
@ 4(ka)* 
= 1Sikd—12(kd)? 
32(ka)* 
where ©, does not depend upon d. For this expres- 
sion to yield a good approximation it is necessary 
that ka>>1 and that ka? >>d (these inequalities 
are equivalent to those given in section 1 in the 
description of Fresnel-region operation). To the 
second order in ka, 


arg ®(d) ~2kd{1— (2ka)~?*]. 


Thus the ‘diffraction correction” to the phase can 
be expressed in terms of a small increase in effective 
wavelength, which in this approximation and in 
this example is independent of d. 

It is of some interest to evaluate this result for 
values of k and a that might be considered typical 


10 


of optical cases, even though the formula does not 

apply, or at least does not apply directly, to any 

optical instrument. Lf one takes \=5000 A and a=} 
em, the fractional increase in effective wavelength 

| given by the formula is approximately 6>< 107%, 

— C3) TE, aperture-distribution. A somewhat more 
realistic—and much more intractable “example of 
Fresnel region behavior is afforded by the assump- 
tion of a TE,o-mode distribution in a square aper- 
ture in a conducting screen. This example is sug- 
gested by an experimental arrangement in whieh 
a square horn-lens radiator is fed by a taper from 
rectangular waveguide supporting the TE, mode {8}, 
If the aperture is bounded by |z|=a/2 and |y|=a?2 
in the plane z=0, we may take 


E, 


cos(mrr/a)e, 
in the aperture and E,=0 elsewhere in the plane, 
From (22) and (23) it follows that 


Syw(1,K)=k,f(K) AK, 


Sio(2,K)=h,/(K) A, 


where 


— .cos (k,a/2) sin (kya/2) 
f(K) —¢ ——<° rs 
. (7 /a)- | ke, 

(Here and subsequently unimportant constant factors 
are denoted (, (’, ete., without explicit definition 
in each case.) Hence (32) becomes 


"(WE ro, ¥ po \1F(K)/? . 2 
&(d) = C” (= k2 + es ) : Rk? | eoedkK. (56) 
} Y Wi ° 


This integral has been evaluated numerically. The 
numerical analysis and the programming required 
for this difficult task were performed by Paul F, 
Wacker and William W. Longley, Jr. Some of the 
results are shown in table 1, wherein the quantity 
labeled Ad is calculated in accordance with the 
definition 


Ad = (2k) [arg?(d) —arg&(0)|—d., 
TABLE | 
l d Ad 
em cm m Microns 
0. 6278 60 2 nM. 06 
. 6278 ht) 10 103. S80 
. 6278 30 2 155. 43 
6278 30 10 503. 28 
l 60 2 2. 32 
1 Ho 10 7.5% 
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The fact that this quantity is negative corresponds 
to the general fact that the effective wavelength 
tends to be greater than the free-space wavelength, 
so that the observed phase increases with d more 
slowly than 2kd. 

D. Fabry-Pe rot with line source. In this example 
the radiated field of the radiating system is taken to 
be identical to that of a simple line current coinciding 
with the y-axis, and it is assumed that the receiving 
system is the same as the radiating system. These 
assumptions vield a two-dimensional, pure TE prob- 
lem in which the free-space radiation pattern has no 
Fresnel region. Although only the most readily 
obtainable results are given here, the example already 
affords an interesting illustration of the behavior of a 
Fabry-Perot instrument in the presence of a continu- 
ous angular spectrum. 

The assumed field being pure TE, we certainly 
have Sijo(l,4,)=0; and since E is asymptotically pro- 
portional to r~'*e"e,, it follows from (25,) that 
we may write Sjo(2,k,)=y"', at least for k2<k’. 
We shall employ this expression for the whole range 
of k,, choosing to justify this analytic continuation 
by the results to be obtained in a moment. Using 
the reciprocity relation (20,) we obtain foe the re- 


ceiving system  S),(2,k,)=y(wu)7'np 'Sio(2,—k-) 
(wun ) Hence the transmission integral (46) be- 
comes 
2 reivty- a 
V(d) —¢ Y _ dk,. (58) 
J-. 1—p’e 
The integral diverges for d=0; this behavior is 


attributable to the implicit assumption of an en- 
forced current in the radiating element (as in the 
dipole example above). If we assume 7 and p con- 
stants independent of k,, d>0, and |p'<1, then (58) 
may be written 


¥(d)—Cr Hdd fe, 


x 2n in 
Goad p / 


ti 


s Lyi (2r 


By introducing a new variable of integration a, such 


that & sin @ y, one may transform the integral 
to one of the standard representations of the Hankel 


function [12] of the first kind and order zero. Thus 
V(d)=C' PD 3p" [(2n + kd). (59) 


nO 


At this point we note that for p=, Wid) becomes 
proportional to //)" (kd), as might be expected for 
free-space transmission between systems of the des- 
cription considered, 

For p approaching unity, features peculiar to the 
Fabry-Perot may be expected to appear—-and the 
series becomes very slowly convergent. A thorough 
discussion of (59) might be premature and will not 
be attempted here (one might wish to consider a 
more realistic model of the Fabry-Perot, for example). 
We shall assume kd>>1 and approximate the Han- 
kel function by means of the first term of its asymp- 
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totic expansion. Instead of (59) we then have 


V(d)~C'' Pp (kd)~*F (p, kd), (60a) 
where 


F(p, kd)=>3 (2n- 1)~3(pet**)2 a 


n=0 


(60b) 


In (60b) a phase angle in p obviously is equivalent 
to an additive constant in d; having noted this, we 
assume the p? is real and positive. 

For p’<1, the series (606) converges for all values 
of d and |F(p,kd)| has maxima for kd=m-z, where m 
is an integer (for p?=1, the series still converges for 
kdAmr but diverges infinitely at the points where 
kd=mr). These maxima are certainly the principal 
maxima and correspond to the passage of axial rays; 
it is not known whether there are subsidiary maxima. 

Since |7\?<1—!pl|* (the equality holding for a 
lossless etalon), + and hence ¥ must approach ze7zo 
as p’—>1 except possibly at the values foe d for 
which the series diverges. An estimate of F(p,mr) 
as p’—>1 indicates that WY must approach zero at 
these points also. This decrease of transmission 
at the maxima as p*—>1 differs from the result 
given by the optical formula (45b) and may be 
surprising. It can be explained in terms of increasing 
selectivity for axial rays, such that a decreasing 
portion of the incident spectrum, and hence de- 
creasing power, is transmitted. 

Finally it may be observed that F(p,kd) is a 
periodic function of d, so that the factor d~} in (60a) 
gives the general trend of ¥ with d. In the optical 
case, WV itself would be periodic. 


4. Appendix 


A. Reciprocity relations. Let E’, H’, and E’’, H’”’ 
denote any two electromagnetic fields (of the same 
frequency) that can exist in the interior of the 
waveguide-space transducer considered. In_ virtue 
of the hypotheses imposed in section 2, the Lorentz 
relation ; 


E’’< H’)=0 


v-(E’ x H”’ 
will hold in V{6). 


theorem, one has 


Therefore, using the divergence 


(EX H’’—E’’ XH’) -ndS=0, 


where nis the inward unit normal on the boundary 
of Vand the parts So, Sj, etc., making up the bound- 
ary of V are as defined in section 2. Now, the contri- 
bution of the integral over Sj vanishes, since the 
integrand vanishes there. By using the asymp- 
totic form of the fields for large r, it can be shown 
that the S.-integral vanishes for fields generated 
by any distribution of sources confined to a sphere 
of finite radius centered at O. This restriction 








apparently would rule out excitation of the system 
by plane waves incident on S,, which we wish to 
consider, but since we can imagine a plane wave 
approximated arbitrarily closely by a source of 
finite dimensions at a sufficiently large distance, 
there is no real limitation. Hence we may employ 
the relation 


| (E’ « H’’—E” «H’) -ndS=0 (61) 
Sot, 


without explicit restrictions. If in this expression 
one replaces E’, H’ and E’’, H’’ by their representa- 
tions on Sp) and S,, equations (1) and (12), one finds 
after some analysis 


no(atbt’—a'b)) — i) ¥7,,(K)[a’,(K)b!(—K) 


—a;, (K)b},(—K)|\dK=0. (62) 
In obtaining this result the use of the integral 
representation (17) of the 6-function is helpful. 

We now assume that E’, H’ and E’’, H’’ are the 
fields corresponding to excitation by the following 
particular sets of incident waves 


” 
Qo = 0, 


a,(K)=0; = a, (K) =6,.,6(K—L). 


From the scattering equations (18) we obtain 


b},.(K)=Syo(m,K), 66 Sy, (n,L). 


Upon substituting all these quantities into (62) and 
observing that 7,,(L)=7,(—L), one obtains the 
reciprocity relation (20) of the text. 

In a very similar manner one can derive (21) of 
the text. 

B. Two-dimensional formulation. For the discus- 
sion of two-dimensional cases many, if not most, of 
the formulas of the text require modification, and 
almost all the modifications may readily be obtained 
by reduction from formulas given in the text or by 
parallel development. Certain key formulas and 
minor subtleties are discussed here. 

It is assumed that all quantities are independent 
of the transverse coordinates y and k,, so that these 
coordinates will be absent from all formulas. 

Under the above assumption the unit vectors «;, 
«x, designating the “parallel” and “perpendicular” 
electric-field directions degenerate to «,—e, sgn k,, 
K,=e, sgn k,. The inconvenience of the sign re- 
versal is avoided by adopting e,, e, as the basis 
vectors for the two polarizations. Sums over the 
polarization index m become sums over z- and y- 
components. This change induces a few further 
sign changes, the key one being in the reciprocity 
relation (given below). 








ee 
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The two-dimensional forms of (14) and (15) may 
be written , 


E(r), | Bk, e*+*"dk,» (14,) 
B(k,) = (27) em | Ee het y (15,) 
where r=se,+-2e,, y’=k*—k?, and _ integrations 


with respect to |» and x are understood to be taken 
over the range (—o, ). (The numbering of these 
two equations indicates the scheme to be followed 
in this appendix.) 

Power expressions such as (16) and (4) must be 
reinterpreted as power per unit length in the y-diree- 
tion. Equation (16) becomes 


P= | do 1m| Om (ke) |?dk,- (16) 


It is convenient to match the loss of a factor of 27 
between (16) and (16,) by renormalizing the basis 
fields at Sy so that (4) becomes 
P. TN | 'dy\” by *). (4,) 
With this renormalization the numerical factors in 
the reciprocity relation remain unchanged. 
The reciprocity relation becomes 
NoSo1 (Mm ke) = nm Sio(m,—k,)> (20,) 
This is probably best established by a derivation 
parallel to that used for (20). The disappearance 
of a minus sign between the members of (20.) as 
compared with (20) is due to the adoption of e, 
and e, as basis vectors for the two polarizations. 
In the two-dimensional case the asymptotic form 
of E(r), for large r is 


, Qn } eo : ‘ 
E(r):, asymp -(=) k cos 6’r- 2e"’ Bik sin 6’), (24) 
where @’=tan7!(z/z). From this follows 
. 1 
Ss I uk ¥ 1 1» ikr A »~ 
. 10 z) =) doy re E(k; kh), asymp (=02) 
for k?<k*’, where k=k,e,+-ye.. 
The power radiation pattern becomes 
p= (e/u)*k cos’6’b-b (26,) 


where b, the complete vectorial spectrum in the two- 
dimensional case, is b=6,e,+)b.e. and 6,=—k,b,y"". 

The two-dimensional forms of equations such as 
(18), (82), and (46) may be written without difficulty. 
(The equation corresponding to (32) is written as (53) 
in the text.] 


(1) 


(2 


) may 


(14,) 


(15,) 


itions 
taken 
these 
owed 


st be 
| 1rec- 


(16,) 


if Qe 


basis 


20,) 
tion 
nee 
) as 


ee, 


rm 


4.) 


1] For 
[1] 7 


(2) 
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Some Solutions for Electromagnetic Problems 
Involving Spheroidal, Spherical, 
and Cylindrical Bodies 


James R. Wait 


(September 29, 1959) 


Solutions are presented for the low-frequency electromagnetic response to an oscillating 
magnetic dipole by conducting bodies of simple shape. The quasi-stationary approximation 
is employed throughout, which is valid when the relevant dimensions of the problem are 
all small compared to the free-space wavelength. This amounts to matching solutions of 
the wave equation within the bodies to solutions of Laplace’s equation outside. The results 
have application to geophysical prospecting. 


1. Introduction 


Electromagnetic methods of geophysical exploration utilize the fact that the conductivity 
of massive ore bodies is much greater than the surrounding barren rock. The general scheme 
is to set up a primary or exciting field by a current-carrying loop and then to detect the second- 
ary field or response of the body by means of a receiving loop. The operating frequency 
should be sufficiently low that the attenuation by the surrounding barren rock is negligible. 
This usually requires frequencies in the audio range. The literature on the subject is extensive 
and here only certain representative papers are referenced [1—4].' 

Theoretical approaches to the subject are usually restricted to highly simplified situations. 
For example the exciting field is often assumed to be uniform or the body is taken to be perfectly 
conducting. While such solutions do indeed provide much useful information, there is a need 
to consider situations of a more general nature. For example, in the practical methods of 
electromagnetic prospecting the transmitting and receiving coils or loops may be located at 
distances from an ore body which are comparable to its maximum dimension and thus the 
uniform-exciting-field assumption is not valid. Furthermore, the frequency which is of the 
order of 500 eps is sufficiently low that the ore bodies seldom behave as if they were of perfect 
conductivity. For the above reasons, it seems worthwhile to set up solutions for certain 
idealized cases which do not suffer from such over-simplifying assumptions which are usually 
present. The geometrical forms considered are the prolate spheroid, the sphere, and the 
circular cylinder. Because of complexity the spheroid is taken to be perfectly conducting. 
The sphere and cylinder are assigned a conductivity and magnetic permeability which are 
finite. In each case the solution is presented for the case of an arbitrarily located magnetic 
dipole. The results are in a form which is suitable for computation. 

In each case treated, special attention is paid to the equatorial plane which contains the 
source and observer. In the three instances, the source is at C and the observer at 7? as 
indicated in figure 1. For a (y-directed) magnetic dipole of strength K at C, the primary 
field at P is 


Ad 
Hi=—>; : ; (amp/m,) 
‘ “Ts 


where Kadi is the magnetic moment. Note that Kd/=(amp-turns) < (coil area). In what 
follows the secondary fields are expressed in terms of x (distance from coil axes to center of 
body), y and y’ (coordinates of source and observer), s (distance between C and P), ¢ con- 


ductivity of body (mhos/m), and u/yo (magnetic permeability ratio). 


eS 


1 Figures in brackets indicate the literature references at the end of this paper. 
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~<t — 
o>} 
Ho (o = 9) 


Figure 1.—Cross section of spheroid, cylinder, or sphere in equatorial plane (z=0), showing transmitting and 


receiver coils. 


2. Low-Frequency Electromagnetic Response of a Highly-Conducting Spheroid 


Prolate spheroidal coordinates are introduced (p, , 6). They are defined by 


9 9 
. ~< 





=1 (eq of spheroids), 





e(m?—1) "em 


(1) 


9 9 
« ~~ 
~ 


p 
se t=! (eq of hyperboloids) 
eI) Tee ioe ae 
where p and z are the usual radial and axial coordinates in a cylindrical coordinate system. 
It thus follows that 
p=c|(1—6?) (n?—1)]} and z=cn6. (2) 


where ¢ is the semifocal distance. These prolate spheroidal coordinates are taken to be con- 
focal with the spheroidal body whose surface is defined by »= 1 (see fig. 2). 


b=+1 major axis = cn 


o 
conducting spheroid n=, | 2 minor axes = a=c J ” -1 
oO 





oe 
c 
t 

















6=0 P 
n=const. 
6 =const, 


Figure 2.—-The spheroid and spheroidal coordinate system 


In the region exterior to the body (i.e., 7 >) the fields are a solution of Laplace’s equation 
if the frequency is sufficiently low (i.e., significant distances should be much less than the 
wavelength). Furthermore, in this case the fields can be derived from a magnetic potential 
2. Thus, 

H=—grad Q 
where 
v72=0 for ">No. 


Now, solutions of Laplace’s equation in spheroidal coordinates are well known and are of the 
form (for integral values of m and n) 


[AP% (6) + BQ? (6)| [A’P2 (n) +- B’ Qi? (n)] e*'™*, (3) 


where A, B, A’, and B’ are constants and, following Smythe [5], P™(u) and Q"(u) are associated 
Legendre functions of argument 4. They are defined by 





ra (us) 7 (u?— 1 ym od / n (H) 
du 
for u>1 4) 
Qe (u) = at@—1) 2 Ga) 
du 
or 
mD /..\ 1 
Pr. (qm) (1—p2ymr2@ I 2) 
du 
lintel (5) 


Qn (u) (1—p2ymi2 27 Qn (H) 
du ; 





in terms of the (ordinary) Legendre functions P,(u) and @Q,(u). It should be noted Hobson 
introduces a factor (—1)" on the right-hand side of the latter two equations and his definitions 
are thus slightly different. - 

At the point C it is assumed that there is a magnetic point charge of strength, K (see fig: 
3). The magnetic potential at P can now be written 


a 
I= 7 pte . (6) 


Ci(n', 6', o') 
R 
P(n, 5, 9) 








No 


Ficure 3.—Conducting spheroid and source C and observer P. 
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where ? is the linear distance between C at (n’, 6’, ¢’) and P at (n, 6, ¢) and Q* is the secondary 
influence of the body. To express 2 in prolate-spheroidal coordinates, we make use of the 
inverse-distance formula 


. Bee. Mh. ; (n—m)VFP i 
>==- >> D> (2—6,)(—1)™ p | (22+ DD Pra )Pr(n)P2(6’) cos m(e—¢’), (7) 
R cx m=0 (n+ mM)? 
which is valid for 7<7’.[6)=1, 6,=—0 (m#0)]. 
The potential Q is a solution of Laplace’s equation and since 1/2 is also, it follows that 
V?Q'=0 for 7 >. Thus 2* must contain terms of the type 


(2? (9) Pe (6) cos m(o—¢’) 


only since P'(n) is infinite at yo and Q(6) is infinite at 6=+1. The potential exterior 


to the spheroid is thus written, for »<n’, 


K 2 a ie 
=) 1S) [MnP (n) + Bran Qe (n)|P2 (8) cos m(¢—¢’), (8) 
4 n=0 m=0 
where 
l (u—m)! , 
Mian=— (2—8n) (—1)"(2n+ 1) | HI On (n’) Pr (6"). (9) 
c (im): 
The unknown coefficient is now found from the boundary condition, that on the surface 
of the spheroid the normal component of the total magnetic field must vanish. Thus, 


022 
—=() at >No. (9) 
on n= 
This is satisfied if 
Pr’ (no) 
... M " (10) 


mn Qn’ (no) , 


where the prime indicates the derivative with respect to the argument of the Legendre func- 
tion. The solution is thus given by 


oar (no 
. i> >> Y Mes mi? a y &n (n) P" (6) cos m(o—¢’) (11) 
R n=0 m=0 (J" 


valid for 7 >. To simplify things a bit we will consider our source to be a y-directed magnetic 
dipole T located in the equatorial plane of the spheroid (z—0) (see fig. 4). 
The potential at P? (x, y, z) is thus given by 


O 
a: a 
dl oy 9 


n=0 m=0 


Kal Q wt, (nd s 
=+ Ag {- Oy’ pt : ee mi! Qn (n) 


OM» 


dy’ cos m(@ saik cos m(¢- 6) Moun | ° (12) 


Pro) [ 


The magnetic field at 2 is obtained from 


H crad ®. 


g 
In the case of the y-component, for example, we have 


OP orn - 
a> = dy’dy (13) 
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(X>y, z) 


C (x', y', 0) 





iy 








FiGuRE 4,.—The conducting spheroid with locations of source C and observer P. 


But 


2 _I2 dn, 92 da» 


dy On dy Od Oy (14) 


and therefore 


ero a2 2s) $t) on (dO¢ Oo") 4 02 Sry 3¢ ;) o’?2? On’\ P 
Sydy" Oe’ Oe (By (5 t55'd6 5a) 5 )t 0¢’dn \dy/ \Oy’ t5," 0d (5%) 5) (15) 


Since 





y=cy (1—6*) (n?—1) sing 


and 





r=cy (1—6") (n?—1) cos ¢, 
it readily follows that for 6=0 (or z=0), 


On_sin¢ yn— l and = et Sm 





——— 16) 
oye oy PI : 


and similarly for the primed quantities. Thus, 


02 
H,=dl dy’dy 
for a y-directed dipole in the plane z=0, 
K 
Q= 40° 
inh . 








where 


K 2 (7—@ Pr’ (no) 
= — nm(—1)"(2 Dm ( s m m 
on 24 p> Em ( ™(2n rye ; iE Qn’ (no) [P™(0)|?cos m(@—o Or (n OR (n), 
where 
Em = 2—6,,=1 for m=0 
2 for m0 
and 
6=6’=0. 
Now 


o* , 1 , rm , , 
dydy? 20 Qn ) cos m(¢—¢’) =V2(n,n’ ,o—¢’) 
*—1)\ (sin ¢’ ¥(n’)?— 1 
n )( c n’ ) 


+ Q(n)Q(n’) m? cos m‘¢- #1) (8. : -)(# 4 -) 


v7 © y(n’)? 


—Q(n)Q’ (n’) m sin — 4 \(AF VG! i ') 


c 5 


—Q’ (n)Q’(n’) cos m(¢- o) (os 


(17) 


+-Q’(n)Q(n’)m sin m(o— oy 22 ons as -) (18) 
, “—] 


. a Cc v(m’)? 


where Q=Q”". 


Finally, 
H,=H? + H; 


where 
yp — Kal re 1 _wdl 3(y ~y" )? he 4 
4 oyoy’ R 4 PS 23 
Kal . 
ea for 2. 
2rs 
Noting that 
R= y’ —y| == y—y’ 
and 
K : | (n— m)! card (no) - 
H; — : > — ad 2 - m m 
’ dnc 2, em 1) ( “ rl) (n rs ‘m)! og ( No) -(P2 (0) )}? V3 


where V® is as given above. 


For the special case (i.¢.,2=2’ =0,|y’ —y|=s, and r=z’) 
as shown in figure 5, we have 


a@=c(ni—1), a=eyn?—1, 
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(19) 


(20) 


y (x,y,0) or (n, 0,6) 
Field Point P 


s (y'-y=s) 


Source C 
(x, y',o) or (n', 0, o') 














Ficure 5.—Special configuration of (dipole) source C and observer P when they are located in equatorial plane of 


spheroid. 


where @ is the semiminor axis, ¢ is the semifocal distance, 


pP=C(n*—1), p=eyn*—1, 
and 
(p’)*=e"[(n’)?—1], p’=¢y (n’)?—1. 
Also it should be noted that 
s=[p’+ (p’)*—2pp’ cos (¢—¢’) |}, 


y’ D 
tan ¢’ y and tan @ . 
I Jt 
If } is the semimajor axis, 


b=cn. 


3. Low-Frequency Electromagnetic Response of a Finitely-Conducting Sphere 
With Arbitrary Permeability 


Spherical coordinates (r, 6, @) are introduced as indicated in figure 6. The surface of a 
spherical body (of conductivity o¢ and permeability ») is defined by r=a. Exterior to the 
sphere the fields can be derived from a magnetic potential Q. Thus, 


where HT: —grad Q 
for r>a. 


vV72—0 


Inside the sphere the fields are a solution of the wave equation and are best derived in terms 
of scalar stream functions (following Schelkunoff) [6]. In the present problem we need only 
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6=0 P (r, 6,0) 


fo 


, c(r', 8, ¢') 








x 


FIGURE 6. The sphe rical body and source C and observer P with the sphe rical coordinate system. 


consider the TE modes which are associated with the induced magnetic dipoles. Thus, 





5 
l oO ' 
H=~~ (53-7) ¥ 
l ere | 1 ©& 
iuwr E 6 O¢ (sin 6 oO” t+ ] 3 | v 
for r<a, (21) 
1 Ow 
IT, Lon Brae 
l oy 
H, imwr Sin 6 Orde , 
where y is the stream function and 
=| tpw(o+ we) t~ (iopw)?s (22 
whereas 
es 02 7 
H~,=— = »=-— 
Lp or v or 


I a tare 1 2& 
Legal | a 6 O08 | mv? O60 + 6 =| Yo 





>forr>a, (23) 
1 OW 02 
He iuwor Or00 roo 
l OY ov 
H,=- = — , 
iywr SIN # Ordd r sin #0 5 


where Yo is the stream function pertinent to the exterior region and is related to the magnetic 
potential by 


o) OY 
ium OF 


22 


(24) 


The stream functions satisfy 


oo |. 45 o l y=7 ry forr <a 
ee 4 e Or 
: IT ae 6 O86 ( ain @ so) tanta & 0 O¢” w—0 for r>a. ana 
Thus they are of the form 
I, 
- Ww) P™ (cos @) for roa (26) 
K, (yr) 
and 
e n+1) Px (cos 8) for r>da. (27) 


a a. 
In the above /, and A, are the spherical modified Bessel functions defined by 


rn. (—1)"(n4+m)! _— n+m)! 
[ p aE aes +(e J —s (28) 


at m!(n—m)'(22)™ »m'!(n—m)!(22)™ 


vA 


I,,(z) 


il 


and 


" = (n+-m)! 
K,(z)=e—? >* a 
7 my m'(n m)!(22)™ ) 
following Schelkunoff {6}. 
Now for r>d. 


Q=2?+-0%, (30) 


where the primary source is denoted Q? and the influence of the sphere by Q. Assuming a 
magnetic charge (strength K) - ated at point C with coordinates (r’, 6’, ¢’), the potential Q? 
at point ? with coordinates (r, 6, @) is given by 


i! 
7 4rh 
K = (ni m)! - 
oe 25 Dy Em 7 rt (p*yati /’m (cos 8)Px (cos 6’) Xcos m(¢—¢") (31) 
47 n=0 m=0 nam)! (rv 


for r <r’, using a well-known addition theorem for the inverse distance in spherical harmonics 
[5|. This suggests writing 


K =~ (n—m)! ie : q?**+! - . : 
" dar D> 2 Cm (hi Mm \! ( n ] ) Sn (7 ut 1pm S17 rs (COs 6) P* (cos 6’) cos m (o—¢’ ), (32) 
. n=0 m=0 T . 


where S, is an undetermined coefficient. The appropriate form for Yo is thus 


l 


Kiu Ma. ah (n—m)! r\er q?*t! : Pr (cos ona (cos 6 
to ——— 2, 2 Oo (37) — Gyre Ss 2 Xeos m(6—¢"). 


dr n=0 m=0 (li- m)! +1) 
(33) 
For the interior the solution must be of the form 
PW 9 n ) a” a” n+1 
Ap owe (n—m)! 7, (yr P" (cos ¢)P* (cos ¢’ ) B.C: ; (34) 
¥ et es Em m , 
4 n=0 m=0 (n- m)'T,, (ya) (n+-1) 


A 
where the Av, (yr) solution has been rejected since it becomes infinite at r=0 and where B, is 
an undetermined coefficient. 


23 








The boundary conditions require that the normal flux (u//,) inside the sphere is contin- 
uous with the exterior normal flux (u//,) at r=a. Also the tangential fields 7/7, and H, 
must be continuous at r=a. These conditions are satisfied if 


oY] wia™ 
*s. =o or (0) 

and 
v=: : (36) 


This leads to 


4A, (a) — (n+ 1) (u/p) 


S.=— - (27) 
A,, (a) 4+ N (L/L) se 
where 
dl (a) 
ee 
da d 5 
A, (a) =—— =a -- [log, J,,(a)] (38) 
I, (a) le : 
with a=ya. Also, 
(M/ Mo) 
B,= . (39) 


1 p/ o) +A, (ex) 
The above expressions can now be used to obtain an expression for the magnetic potential 


due to a dipole for r>a in terms of the potential Q due to a point charge. For example in 
the case of a y directed dipole at (r’, ¢’, 2’) we have 
, ie 
P, —<dl ; Q)- (40) 
oy 


The corresponding magnetic fields are thus given by 


> 
H=—grad (41) 
and the y-component, for example, is 
02 
H,=dl —* (42) 
y oyoy’ / 
Again we shall consider the equatorial plane (see fig. 7) which is z= 2’ =0 or 
6=0’=7/2. Here r=r cos z’=r’ cos ¢’ 


y=r sin @ y’=r’ sin ¢’, 
and thus, 


syasin 6° +e0s 6.53 (43) 
and 
o 5 oe 
aye dp? t 008 5 Oe” (44) 


The field is thus obtained conveniently from 


[si a ey ae o> — >? 
H, dl sing sin o> 7op tT sin & mes +cos @ aeha!”, Top t OOS ated ; |2. 
(45) 
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H, 


6) 


‘) 


It should be noted that 





g’ =r’ for 6’ = 2/2 
(46) 
p=r for 6= 2/2 
Writing 
H,=H?+ Hi}, 
it is seen that 
HT? = | a a where R=¥V(r—2’)?+ (y—y’)? for g=-2'=0 (47) 
4m Oydy’ R i ’ 
and thus, 
Kdl 0 y-y" +o )? K 
H?—~——_ — (> )=—dl| — 
a na le Rs at R* 4x 
Kal ? 
— = where s=|y—y’|. (48) 
4rs vil 0. 
The secondary field is derived from 
Ke 2 (n—m)! a’"™*'n , Pt (cos 0)P% (cos 6’) 
05: Fs y ” 7 Dy . s — i 
de <—, “™ (n t m)! (7’ persgers ; n+1 — m(o %) (49) 
Finally (for r2 a) 
Kw 2 (n—m)! es” 
ae SS - §,{[P™(0)2U wm —peccy 5 
Hy, rap >p 2e ™(n nen Org nln) PU, (7? )etipeti (50) 
where 
1 
ins = l(n +1)? sin @’ sin @ cos m(o—¢’) 
+ m(n+1) sin ¢’ cos @ sin m(@—¢’) 
—m(n-+-1) cos ¢’ sin @ sin m(¢—¢’) 
m’® cos ¢’ cos ¢ cos m(@—¢’)), (51) 
where 


tan = LY, 
ya 
, 
tang’ —". 
- 

The quantity S, plays rather an essential role in the final result. When the conductivity 
is great or at high frequency, the argument ya of the Bessel functions may be large. In fact 
if \ya| > >1, S, approaches unity for alln. The behavior of S,, as a function of ya and y/yo, is 
a convenient way to illustrate the frequency dependence of the secondary fields. Numerical 
values of the real and imaginary parts of S, are given in some detail for n=1 in previous 
papers [2, 8] and less adequately [9] for n=2, 3, and 4. Fortunately, if the source C or the 
observer P (for any respective orientations) are a distance from the center of the sphere which 
is large compared to the radius, it is seen that only the term corresponding to S, is important. 
The response of the sphere is thus proportional to the quantity 


owt. 
S,= yal; (ya) — 2KT, (ya) where K=p/Ho. 


yal (ya) +KT,(ya) 
The real and imaginary parts of S, are plotted in figure 8 as a function of the “relative radius” X 
defined by 
=ViX or X= (ouw)*a. 
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For small values of Y corresponding to the lower frequencies, S, is negative and, in fact, 
it approaches the value —2(K—1)/(K-+2). This could have been derived directly from 
magnetostatics. At the higher frequencies S,; approaches unity which corresponds to no 
penetration of magnetic flux into the body. 


4. Low-Frequency Response of a Conducting Cylinder With Arbitrary 
Permeability 


We now wish to calculate the electromagnetic response of an infinitely long conducting 
cylinder with any permeability. (See fig.9.) Exterior to the evlinder p >a, the fields are a 
solution of Laplace’s equation and are derivable from a potential. Inside the cylinder the 
fields are a solution of the wave equation and are derivable from a magnetic Hertz vector II* and 


an electric Hertz vector IT with only z-components. Thus (for r <a) 


H—=(—y*+ grad div) lI*+—— eurl TI, (52) 
Uw 
where 
I1*— (0,0, 11") (53) 
and 
II— (0,0,11.). (54) 
Solutions of the wave equation 
» » 1 * _— 
(v?—+7’),,7 =0 (99) 
are of the form 
I, (ap) 
ap ime, —ith (56) 


K,, (ap) 


where a=(y?+h?)? and J,, and K,, are evlindrical Bessel functions of the modified kind defined 


by G. N. Watson [7]. 











FIGure 7.—Special configuration of (dipole) source C and observer P when they are located in equatorial plane of 


sphere. 
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Figure 8.—The variation of the response of a conducting sphere as a function of X(=(oepw)*a) for various 
permeability ratios K(= p/p). 


The real and imaginary parts of S; are shown. 
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Figure 9.—The cylindrical body and source C and observer P with the cylindrical coordinate system. 


External to the evlinder (7 >a) 


H=—grad 2 (57 
where 
S22" -+- i. (58) 
For a magnetic charge at €' (p’, o’, 2’) and the observer at 2? (p, @, 2), 
Pee (59) 





P(x, y,0) 


®\ PHC (x, y', 0) 
GPa 











Figure 10.—Special configuration of (dipole) source C and observer P when they are located in equatorial plane 
of cylinder. 


where ?? is the usual linear distance between C and ?. Using a formula given by Watson [7] 
we find 


eo go ots , 
>>=- | Ky (hp)e 2 dh, (60) 
 @2c. 
where 
R=, (p)?+ (2z—2’)? (61) 
with 
p=|p’+(p’)*—2pp’ cos (¢ ¢’)}?. (62) 
Using the addition theorem 
Ky(hp) = >> é€nin(hp’)T, (hp) cos m(o—¢’) (63) 
m= 
for p<p’ and with 
@=1, @n=2 (m=1,2,3 ... ), (64) 
it now readily follows that 
x 2 fy. , hie—ar im(o~<¢ 65 
0° = ph K,, (ho’) { Im(hp)e~ 2? dh} e— im e-©, (09) 
4 M=—@ .J —@ 
This latter equation can be written symbolically in the manner 
0?—TT,,(hp) (66) 
where T is the operator 
x @ ft... pie ) - 
a3 Ps . mee vi. wee er? Re MCC. (67) 
This suggests writing 
. O=—TA,, (bh) K,, (hp) (68) 
or 
Q=T UT n(hp)+ An(h)Kn(hp)| (69) 
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for the external fields and 
li*=Ta,,(h)Im(ap), (70) 


11,=Tb,,(h)Im(ap) (71) 


for the internal fields. The corresponding magnetic field components are 
H,=—Th[In(hp) +An(h) Kn (hp)] ) 


um 


H, -T = Um(hp) +Am(h) Kn(ho] > rDa (72) 





H,=Tih{In(hp)+Am(h)Kn(hp)} J 


and 
sy 


H,=T {. thoa,,(h)I’.(ap) —~— ba(h) Ia(ar) } 
pw p 


H,=T { din(h) Im (ap) —2— bn (h) In(ap) } r rsa (73) 
p 1p 





H,=T{ —a2,(h) In (cep) } ) 


4 


The unknown coefficients a,,, 5,,, and A, which are functions of h, are determined from 
the boundary conditions. These are the normal flux, (i.e., w//,) and the tangential fields, 
(i.e., 7, and /7,) are continuous atr=a. The explicit result for A,,(h) is found to be 


~ ~ my" l 
I, (ha) — (p/m) | Fa(aa)+ mit ~ ] 
awd J,(aa) I, (ha) 


A,, (a4) =—} - 7 , (74) 
= ~ a ~ m*y* l K,,, (ha 
| Katha (u/ ty) | Fa(aa)+ L | 


ij,2.4 > 
a‘h'a I, (aa) 


W here 
~ T’,(a) 
I, (a) se 75) 
al, (a) ( : 
and 
a, K°, (a) . 
K,, (a) = (76) 


= > . 
ak m(Q) 


The preceding results can now be employed to calculate the magnetic field of a magnetic 
dipole. In the case of a y-directed dipole at C’ the potential at Pis —dloQ/dy’. The y-com- 
ponent of the field at ? is obtained from 


ee is. 77 
u i Oydy” (4 i) 
Restricting our attention again to the plane z=2’=0 and for z=2’ (see fig. 10), we have 
bull ” 
a = (78) 
where s=y—y’ >0. Noting that 
, o o 
sin @ cos } ~y (79) 
—? poo ' 


it readily follows that 


HT’ > > 


Kdl ts (te m im a> 
-> A,,(h) | hk), (hp) sin @— wn K,, (hp) cos @ 
J-4 p 


4r* m 


| AK (he) sin 745 K,,(hp’) cos |. MO~- dig ™O-e), (80) 
p 


where A,,(h) is given explicitly above. 
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5. A Related Cylindrical Problem 


The above general solution for an infinitely long cylinder with dipolar excitation is rather 
complicated. The presence of the infinite integrals and the existence of Bessel functions of 
generally complex argument complicate the numerical aspects of the problem. It appears, 
however, that simplifications can be made under certain limiting conditions. For example, if 
the radius a is small and at low frequencies, ha is much less than unity over the important range 


of the integration. Thus, 


, mKT,,(ya)—(ya) I}, (ya) I, (ha) 
An(h) ~—| \ 2 = ] - (81) 


mKT,,(ya) + (ya)I/,.(ya)_| K,,(ha)’ 


where K=p/y. In the case of very high conductivity |ya. -1 this simplifies further to 


T,,, (ha) 


wi , 9 
K,, (ha) (82) 


A, (h) ~ 

The square bracket term in equation (81) thus characterizes the frequency dependence of the 
response and is a function only of A and ya. 

The secondary magnetic field 77° for the case of finite ya can then be written in terms of 


~ 
value of //* for infinite ya in the following way (for a given mode number m) 


Poses 
—— ol. (m=123,...), (83) 
[HT*| ya->@ 


where 


mKT, (yi X)—yi XT), (yi X) 
ad \ \ \ (84) 


mK, (yi NX) +y7 NI}, (yi X) 


» 1 
and where .Y= (ow) 2a. 


The total field is, of course, a sum over all modes, but because of the factor J,,(ha)/K,, (ha), 
the relative response decreases rapidly with increasing value of m. This follows from the 
approximation 

T,,, (ha) (ha)?” * 
~ (So) 


K,, (ha)~ (m!) (m—1)!2?""! 


which is valid for (ha)’<1. 


To shed further light on the above approximate form, it is constructive to examine a 
related problem. If the source is an infinite line source of electric current J, located at (pp, dp) 
and parallel! to the z-axis, it is known that the secondary fields due to the presence of the evlinder 
are [10] 


I ae a-” 


IT; —<- DS Dn ax SiN M (P6—>p) (86) 
<7 l Py Pp 
and 
, I, si es ies 
HH, — . # maj, €OS MW (o- dy). (87) 


2r m=1 


" pip 
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The structure of this exact two-dimensional solution is very similar to the approximated form 
of the three-dimensional solution mentioned above. Apparently if @ is reasonably small, the 
induced currents flow mainly in the axial direction even for a localized excitation and thus the 
frequency dependence of the induced eddy currents are adequately described by the form of the 
solution of the corresponding two-dimensional problem, 

The numerical values of the real and imaginary parts of the function 7, are illustrated in 
figure 11, where they are plotted as a function of YX, the ‘relative radius” of the cylinder, for 
various values of the permeability ratio A. It is noted that the factor 7,, corresponding to 
the dominant mode, is very similar to the factor S, for the sphere shown in figure 8. 
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Figure 11.—The variation of the response of a conducting cylinder as a function of X(=(ouw)* a) for various 


permeability ratios K(= p/p). 


The real and imaginary parts of 7), 72, and 73 are shown 
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Kantorovich's Inequality ' 


Morris Newman 





30 : 
(August 18, 1959) 
al, ; . ; : 
(An elementary proof with «a generalization of an inequality of Kantorovich is given. 

)): 
in Let A be a hermitian positive definite matrix |  Proor. Since A is hermitian, A is unitarily 

with smallest eigenvalue a and largest eigenvalue 8. | equivalent to a diagonal matrix. We may therefore 

Then Kantorovich’s inequality * states that for all | assume without loss of generality that A=(a,,), 
| vectors z of unit norm, where m<a,<M. Then also F=(f(a,6,) and 
cs | G=(g(a)6)). Suppose that the i” coordinate of zis 

i 1 a\} \3 ) ? x, and set ¢t,—!z,|?, so that 
a. (Ar, r)(Aq'r, n<t{(§) +() } . (1) | , =e 
; a 
n t,>0, 2ts=!1, (6) 
_ In this note we give an elementary proof of (1) ; os . 

which allows an easy generalization. The results | “ being of unit norm. Then K becomes 
. obtained do not depend on the order of A and this . 
’ is accordingly left unspecified. K= Dos (ait Dg @)ti- 

We assume only that A is hermitian, and that the . 7 
eigenvalues of A are contained in the closed interval Let ¢ >0 be arbitrary and rewrite K as 
m<t<M. | 
Let f(t), g(t) be real functions such that K=Dref(a) to - gay) te. 
0< f(t), g()<@ for m<t<M, (2) 
Then 
f()), g(t)” are convex for m<t<M. (3) 


, a 7 , 
2K?=2(Sref(a)t,)? (> = g(a,)t: ) 
Then (2) implies that the matrices F'=f/(A), | 
G=g(A) are well-defined and are hermitian positive | — | 
definite. We shall prove: <Soef(a)t:+>5 = I(ai)te 


THeoremM. Let x be any vector of unit norm. Put 


K=(Fryr,x)(G@r.r). < max (cf(a,) + g(a,) ), 
E c } 
by (6). 


The Nt for every positive é., 
. We now consider the function 
2h? < mar(ef (m)+—g(m),ef( MM) +-9g(M)). (4) | 
ce ; e 
l 
7 = u=cf(t)+—g9(t) m<t<M, 
If in addition f(M)—f(m) and g(\1)—g(m) are of © 


opposite sign, then j 
Because of assumptions (2) and (3) and the fact that 


_¢ is positive, we have that w is convex and positive 


i . l a ar : 
2h? <rf(m - g(m), (9) | for m<t<M/. This implies that 
where max u=max (u(m),u(M]M)) 
m<t<M 


1 

( g(M)—g(m)\2 y 
ahd tent ; which proves (4). 

| If in addition f(M)—f(m) and g(.M)—g(m) are of 
| opposite sign, the choice c=r is permissible and makes 
. The preparation of this paper was supported (in part) by the Office of Naval 1 1 
esearch 
‘L. V. Kantorovich, Functional analysis and applied mathematics, Uspekh rf(m)+—g(m) =rf(M)-g(. 
Vas et nae santa end speed mathemati, Uap f(m) +—9(m) —rf(M)-9(M), 
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which gives (5). This completes the proof of the 
theorem. 


Kantorovich’s inequality (1) is the case f(t)=f, 
q(t) a | for which r=(mM)~** (which is permissible 


when A is positive definite). 

Ky Fan in a written communication to the author 
points out that the theorem can be generalized as 
follows: 

Let A be a hermitian matrix of order n with all its 
eigenvalues contained in the closed interval a<t<B. 


Let 1, %, . . . , ty be vectors in unitary n-space such 
that 
k 
aol |?=1. 
fT 
Let Sis Fe, ss » Sei be positive CONE x functions for 
acts. 


Furthermore if A is positive definite and a> ¢ 
Kantorovich’s inequality can be replaced by 
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i=1 i=] . 


| aw : ai\? h)2 
| 2 (Ar,,20 {> (A roe) } <3{() +(2) }. 


Atisfying 


fi a 


d a>), 


JOURNAL OF RESEARCH 


Vol. 64B, No. | 


{ the National Bureau of Standards 
January-March 


B. Mathematics and Mathematical Physics 
1960 


A Symmetric Continuous Poker Model 
A. J. Goldman and J. J. Stone ' 


August 


Beginning with von Neumann, 


mathematicians concerned with the 


IS, 


1959 


rational analysis of 


conflict situations have realized that investigation of accurate or simplified versions of com- 
mon card games leads to techniques and insights applicable to situations of military or 


economic interest. 


cated than the original von Neumann game is solved. 
permitted. 


): no raises are 


a low 


an ante of 1 unit (a>b>I1 
strategy, which forbids bluffing on 
the von Neumann model 


1. Introduction 


The technical literature of the mathematical theory 
of games of strategy can be roughly subdivided into 
i) general papers, which deal with conce pts, existence 
theorems and computational methods relevant for 
broad classes of games, and (i7) specific papers, in 
which the solutions of particular games are deter- 
mined. Although many of the studie ‘'s in the second 
category have dealt with problems suggested by eco- 
nomic or military conflicts, e.g. [7, 2]?,a substantial 
fraction have dealt with accurate or simplified ver- 
sions of common card games, the best known example 
Neumann’s analysis of a simple poker 
model ({3], sec. 19.4 to 19.10). Bridge and baccarat 
have been treated [4, 5], but poker has remained the 
favorite topic for these studies [3, 6 to 10], perhaps 
because the poker models investigated have pre- 
sented just the right degree of challenge (“not easy, 
but not impossibly hard’’) to the ingenuity of the 
analyst. Of course the techniques and insights de- 
veloped in these investigations are then applicable to 
situations of greater practical importance. 

The present paper deals with a symmetric poker 
game one stage more complicated than von 
Veumann’s original model. The game is continuous 
in the sense that the difficulties [6] due to the enor- 
mous but finite number of possible hands have been 
avoided by assuming von Neumann did) an 
infinite continuum of possible hands. The (more 
“ener ag asymmetric version has also been solved 
by the authors; the results will appear elsewhere [11]. 


being von 


(as 


2. Description of the Game 


The rules of the game are as follows. The two 
players ? and S first ante 1 unit. They then receive 
independent random numbers (hands) from the in- 
terval [0, 1]. Each player knows his own hand, but 
not that of his opponent. The players act simul- 
taneously, each either dropping, betting an additional 


s, Stanford University 
in brackets indicate the literature references at the end of 


Present addres 
* Italicized figures 
tis paper 


hand. 
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In the present paper, a symmetric poker model one stage more > compli. 


There are two bet levels, b, and 
The game has a dees pt Peer 
is shown to vield 


The 


limiting case b= 1 


b—1 units, or betting an additional a—1 units. If 
both drop, then play ends and no payment is made. 
If both bet 6—1 units or both bet a—1 units, then 
play ends and the player with the higher hand * wins 
the ante and bet of his opponent. If one player 
bets a—1 units and the other bets 6—1, then the 
latter can either fold (in which case he loses his ante 
and previous bet to his opponent) or see * by betting 
an additional a—b units (in which ease the player 
with the higher hand wins the ante and previous 
bet or bets of his opponent); in either case play ends. 
Naturally we assume a>)b>1. 

Each player has four courses of action for any given 
hand. These are (1) to drop, (2) to bet 6-1 units 
with the intention of folding if the opponent bets 
a1 units, (3) to bet 6-1 units with the intention 
of seeing if the opponent bets a-1 units, and (4) 
to bet a-1 units. Figures 1 and 2 show the algebraic 
payoff from player S to player R for any choices 
of courses of action by the players. 


A strategy (for either player) is a quadruple 
r(z)=ird2)} (@=1,2,3,4) of functions defined for 
O<2<1. Since r,(z) represents the probability ° 


that the player will choose the 7th course of action 
when his hand is z, the restrictions 


r,(z) 20, ry (2) +o Zz) +973(2) +14(2) = 1 
are imposed. The “smoothness” conditions on these 
functions are (7) that they be integrable, (77) that 
if r,(2") > 0 then, for any interval J of which 2’ is 


an interior ° point, 


| r,(t)dt > 
J 


and finally * (/77) that 


| (ors(e) ar,(t)|dt=0 on [(b—1)/b, (a 


1 )/a} 


} The possibility of equal hands has zero probability and so can be disregarded, 
‘ This possibility distinguishes our model from the case m=3 of the version 
studied in part II of [4] 
More precisely, ri(2) 
6 The words “interior” 
[0,1] 
’ This very weak but ad hoc condition is precisely what is needed to prove 
b) and (c) of lemma 8. 


is a probability density function. 
and **open”’ are to be understood relative to the interval 








implies 


0 [(b—1)/b, (a—1)/a} 


br.(z)—ar,(z) on 


It is convenient to define, for any r(2), 


ad 
» 


r,(t)dt. 


For any r(z), @(2) will denote the expected payoff 
to player /??, using r(z), against a fixed hand 2 and a 
fixed course of action 7 of his opponent. Thus the 
expected payoff function of the game is given by 


PT 


K(r,r*) Dori*(2)G (2) dz (2.1) 


The explicit formulas * for the @,(z) are readily 
found from figures 1 and 2: 


\ 
. Ss l ») ) 4 Ss ] ») > 4 
- ~ » 
R RN 
1/0 | l l 1 O I | ] 
o 1% h i) b 14 bh b h 
an b b a a l h b a 
4/1 h a a 4/;1 bh a a 
FIGURE | Figure 2 
P has the higher hand. S has the higher hand 
G,(2)> | “[r3(t) try(t)+r,(t) |dt 
J 0 “ 
+ | [ro(t)+-73(t)+-ry(t) \dt, (2.2) 
G,(z) | r, (t) —br.(t) —brs(t) + bry(t) \dt 
J 0 - 
+| [—r,(t) +br.(t) +-brs(t) +br,y(t) \dt, (2.3) 
G,( 2) | —r,(t)—br.(t) —br,(t)—ar,(t) \dt 
J - 
+ | |—r,(t) T br,(t) T br,(t) T ar,(t)\dt, (2.4) 
G,(2) | | —1r,(t) —br.(t) —ars(t)—ar,(t) \dt 
J 0 
"1 
+| [—r,(t)—br,(t) +ary(t)+ary(t)\dt. (2.5) 


* Note that the G;(z)’s are shown by these formulas to be continuous functions 


3. Plan of Attack 5. 
It would not be very instructive simply to state 
the solution of the game and demonstrate its corree. req 
ness. Instead, a slightly polished version of the ope 
process actually used to solve the game will be pre- (ba 
sented. This is done in sections 4 and 5, in whieh on 
it is proved (subject to the pleusible assumption 
stated in sec. 4) that the only possible optimal strat. 
egy is the one obeying the conditions (a) to (d) of 
lemma 8, which deal with the behavior of the strat. 
egy on three intervals [0,q), (g.g’), and (g’,1] which 
might be thought of as the sets of /ow hands, jntey- 
mediate hands, and high hands. It is then relatively 
easy to verify that this particular strategy is indeed 
the unique optimal strategy; this is done in section 6, 
primarily by inspection of the expected payoff fune- 
tion (2.1). In section 7 the limiting case b=] 
dealt with, and this game is found to have a unique 
optimal strategy obtained * by ‘‘passing to the limit” 
in the optimal strategy for b>1. The von Neumann 
poker game of sections 19.4 to 19.10 of [3] is shown 
to arise as a “reduced game” from this limiting case 

In section 4 it is shown that 73(2)=0 on [0,1] fo: 
every optimal r(z). The points g, g’ are located, and 
the unique optimal 7(2) is isolated in section 5. The‘ 
arguments of these sections are essentially all appli- 
cations of the following domination criterion: 

The strategy r(z) is optimal only if the following is 
true: For each 2, 7 for which Giz) (as a function of i 
does not assume its minimum, r;(z)=—0. 

The correctness of this criterion follows from the 
smoothness condition (77) and the fact that r(2) is 
optimal only if 


“a 


Is 


cl 
ta 
th 
Ihe 
Ce 


Kirwr)=—min A(r,r*). 


x 


/ 


4. Elimination of One Course of Action 
re 
The elimination is effected as follows. . 
LemMa 1. For eve ry optimal r(z), 
(0,1). 
Otherwise, for some optimal r(z), R,>0 and 
r,(z)~0 for some .. the least) upper 
bound of the set of all such 2, then by (2.4) and (2.5 


r.(z)=0 of 


Is 


(a—b) R,—2b | r.(t)dt<0, 


G,(z’)— G;(z’) 


so that by continuity and the domination criterion 


r(z)=0 in a neighborhood of 2’, contradicting the 


| definition of 2’. 
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The following plausible assumption (which will 
be verified in sec. 6) is now made: 

Assumption: For every optimal r(z), 
(and thus Ry> 1). 


rgl | )> 0) 


This is not quite the case for the asymmetric version see 11 


5. Derivation of the Necessary Conditions 


tO state Since application of the domination criterion 
COrrect. pequires examining the signs of the pairwise differ- 
Of the onees of the G,(z)’s, the formulas for these differences 


%” pre- (based on (2.2) to (2.5) and the condition rs(2) =0 
Nn Which on [0,1]) will first be given: 
IMDtion 
al strat. - 
Strat Go(z) — Gy (z) = (6—1)—bR,— 26 r(t)dt (5.1) 
» (d) of 2 : 
iC Strat. 
| Which ot 
“ee > ‘ : 
S. inter. G)(z)—G@i(2) =bR,+ah,—1—2b | ry (t)dt 
latively 
; Indeed : *2 ns 
ction §. afl ry(tydt (.0 
ff funce- . 
b ] Is . 
unique (7; (2) (i, / ) (a b)R, 2a r(tydt (5.3) 
» limit” we 
umanh - 
Shown G,(z)—G,(2) 1—bR,+ah,—2a | r(t)dt (5.4) 
if Case 0 
0.1] fo : 
ed, and (| >) (7, ( Dh PR, la b ) R, Dh Is ( / dt 
De The a/ 0 
| appli- ; 
w 2a | ry(t)dt = (5.5) 
Wri in 0 
m of r 
Fe (2) 2bR,+-2 Te 5.6) 
ym the (r,4(2 G( PhP, + 2b | ra(ta (9.6 
r(z) is 
The first four of these continuous functions are 
clearly monotone nonincreasing. The following 
table defines the maximal open intervals in which 
these functions are positive (+), zero (0), and 
negative (—). For each interval in the (+) and (—) 
s columns, the result of applying the domination 
10n criterion is given. Inspection of (5.1) to (5.4) 
reveals conditions on the intervals in the (0) column 
() ' Of course lemma | renders superfluous the appearance of r Oin the table 
on 
0 and 
upper 
9°: a» 
(2.5 Difference (+) 
0, 
G,(z)—G,(z) [Og )? 
ro(z)=0 
terion 
y ; y y ; ‘ 
is th Gr,(2)—G, (2) lO.g*! )? 
; ral z)=G 
1 will 
’ G;(2)— G)(2) lO,g**) 
I)70 r(z)=0 
G,(2)—G,(2) lO,g'!)? 


r.{2)==0 





| 
| 
| 
| 
| 


| 
| 
| 


} 
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which are given below the interval. A question 
mark means that the interval may (so far as is 
known at the moment) be vacuous or degenerate 
to a single point; the condition written below such 
an interval is to apply only if the interval is non- 
degenerate. 


The next three lemmas vield the values of ?;, 2, 


and /?, for optimal r(z) (lemma 1 shows that 2,0). 
LemMA 2. For every optimal r(2): 
(a) (.(0) —G,(0) = (6—1) —bR, > 0. 
(b) G,(0)— G,(0) 1—bR,+ah,<o. 
If (a) were false, then the monotonicity of 
G.(2)—G (2) and the domination criterion would 


imply 7,0, so that (a) would be true. If (b) were 


false, then since Gy(1)—G\(1)—0 one would have 
(7,(2’) —G,(2’)=0 
for some 2 in (0,1); by (5.10) 74(2) <0 on [0,27); so 
by (5.4) 
(7,(0) — G,(0) = G,( 2’) — G,(2’) =0 
and (b) would be true. 
LemMa 3. For every optimal r(z): 
(a) G,(0)— G,(0) = (6—1)—bR, = 0. 
(b) (0) —G,(0) 1—bR,+aR,=—0. 
If Gil)—G,1)>0. then the monotonicity of 
G(2)—G,(z) and the domination criterion would 


imply R,=0, so that the one hand the initial 


| assumption would read 


on 


(b--1)—bR, >0, 
(0) ( ) 
(ge gs)? (1)? na 
re(z2)=0 rj(z)=0 (.4) 
fs I nes : A (5.8) 
i hi : RB Ai | (5.9) 
oS: A | B Ay | (5.10) 


7 








and on the other hand P?,;=1—R, (by lemma 1, 
using 2?,=0). Combining these statements vields 
R,>1/b, which together with (5.4) and R,—0 vields 

G(0) G(0) l 


t ak, 0, 


contradicting (b) of lemma 2. Thus G,(1)— G,(1) <0, 
2 implies that 


which together with (2) of lemma 2 


G,(2’) —G,(2’)=0 
for some 2’ in [0,1]. By (5.7) re(z)—0 on [0,2°); so 
by (5.1), 
G,(0) — G,(0) = G,(2’) — G,(2’) = 0, 


and (a) is proved. If (b) were false, then by (b) of 


lemma 2, the monotonicity of Gy(z)—G,(2), and the 
domination criterion, we would have /;—0, con- 
tradicting (a). 
LeMMA 4. kor eve ry oplimal “ata 
(a) Ry, = (b—1)/b. 
(by) R= (a—b)/bla > b) 
(C) R, 2 (a+b) 
This follows from lemma 3 and the corollary 
R,+-R,+ R,—1 of lemma 1. 
Next the points g, yg’ will be located, and the 


behavior of an optimal r(z) on |0.g), (g.g’) and 


(y’, 1} determined. 


Lemma 5. For every optimal r(: 


(2) q) es q (h l } h. 
(b) r,(z)=1 on |O.9g). 
(c) r,(z)=0 on (y, 1}. 
By (5.5) and lemma 4 
G4(z)—G(z)=2b | r(t)dt—2a] ryt)dt (5.11) 


If @i<g'!, then by (5.7) and (5.10) one has ry(2)= 1 
on (g?, gi), m(z)=1 on [0, g%), and ry(z)=0 on 
[0, g‘}), so that the last equation vields 
G,(y'!) —Gi(g") >0: 

continuity and the domination criterion show that 
r(z)=0 in a neighborhood of g'', but (5.4) and (b) 
of lemma 3 show that this contradicts the definition 
of 7’. An entirely similar argument refutes the pos- 
sibility gii<g?!, and so @i=g'!. By (5.7) and (5.10) 
we have r,(z)=1 on [0, g%) and 7,(2)=0 on (9%, 1], 
so that @!=/?,;; an application of (a) of lemma 4 
completes the proof, 


LeMMaA 6. For ere ry optimal rte): 


, 


a l)a, 


(ab) ye (a 


(b) ry) Lon iy’, 
Note first that 7,(2 

on (g°?, 1] bv (5.9), (it) 7) (2) 

and (ii) gg? because Gy(y) 

the fact (see (5.10)) that ry(2 

that 


Non (y, 1] by lemma 5 
G(qg) >0 bv (5.3) and 
Yon |0.q). It follows 


| r(t)dt+ (1—g?)— Ry? 


(ad 


bh). 


The definition of q . together with (5.3 shows that 


| “recat (a—b) Rya—(a—b) ala+b), (5.42 


and combining the last equations vields the asserted 


value of g*. 


LEMMA 7. For eve ry optimal riz). 


G,(2) (i, ) on ld, q'| 
Since ry(2)—O0 on (yg, Tf. (5.11) and (5.12) imply 
Gly’) Gy(y' PAR, (a A\R, (), 


where the lust equality follows from lemma 3. Let 
>* be a point at which G@,(2)—G,(2) attains its maxi- 
mum \/ on {g,g’|; the last equation shows that 1/>0 


If Af >0, then by continuity and the domination cri- 
terion one has 7,(2)—0 (and thus r(z)—1, since 
r(2)=0 on (g,1]) on some interval (2*,2°); by (5.5 
G,(2)—G@,(z) is strictly increasing on which 
contradicts the definition of 2* unless q’, 


which ease M0 by the last equation. Thus 1/=0 
an entirely similar argument shows that the mini- 
mum value of G,(2)— G.(2) on {q.g’ 0, and so 
the lemma is proved, 


Is wi 


The completion and summary of the preceding 
material is as follows. 


LemMMa 8S. Let q (h—1) 6 and y (a—I1)/a, If 
the assumption at the end of section 4 is correct, then 
the only possible optimal strate qy is the strate ay defined 
by 


(2) r,(2) l on [O,q). 
(b) r(z)—a/(a+b) on (g,g") 
(c) ry(z)—b (a+b) on (a.g’ 

(dl) r{z)=] on (q’ 1). 


| 
| 


From lemma 1, lemma 5, and lemma 6 one has (a), 


(d P and also 


on (g.q’). 
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lon (y"?, 1} since (1) ro(z)=—: 


wl 


thi 


6. 


fir 
H. 
if 

se 


mM 


From lemma 4. (BD), and lemma tone has 


‘[br,(t) ar,(t)\dt 0 on (g.q'), 


which by smoothness condition (777) (see section 2) 
"»(2 oo 
I's ~) () implies 


mma 5 
3) and 
follows 


brs(z) airs(2) 0) on (y.q’); 


this and (5.15 imply (b) and (e 


6. Proof That the Conditions Are Necessary 
and Sufficient for Optimality 


WS that : .; ; , 
In this section the notation 7(2) is reserved for the 
particular strategy described in lemma 8, while arbi- 
(».12 trary strategies are denoted r*(z2). 


LEMMA %. Kaw min Wirwr*) if and only if 


sserted 


y* (> obeys 
) a, () on |g f 
2 * , 
by ry(2 r3(z)=0 on (7d 
c ri ( rs(z () on (9,1). 


imply 
The Kiw*) of (2.1) is first written as the sum of 


these integrals: 


7 Let . 
 mMaXt- | G,(2) +3 (2) (G@ey( 2) — G,(2)) + F(z) (G(2) 
M>0. . . 
ion ¢ri- G,(z)) 4+-r¥(2)(G,(2) —G >)) Md. 
since 
(5 5 ’ ’ 
y _—_ | (7,(2 r(2)(G,(2 (75(2)) +-9®( 2) (G3(2) 
which . 2 
gy’, (7.(2)) +-9¥(2)(G4( 2) —G.(2)) |dz. 
M=0 
/ mini-) 7 
nnd so | (ry(z) 4 re (2) (G,(2 (r4(2))- rs (2) (G@(2) 


, Gi (2)) 4-73 (2) (G4(2) Giy(2))\dz. 
ceding 
lt is easily verified that r(z) has the properties 
a. If deseribed in lemma 4, and since rs(z)—0 on [0,1], 
t. then (2) obevs (5.1) (5.6). Using these facts, one finds 
defined that on [0,g) the coefficients of r¥(2) and r¥(2) in the 
first integral vanish, while that of r*(z2) is positive. 
Henee r*(2) minimizes the first integral if and only 
if it satisfies (a). Similarly, 7%(2) minimizes the 
second integral if and only if it satisfies (b), and 
minimizes the third integral if and only if it satisfies 
x 


THEOREM. 7 


is the unique optimal strategy. 
“us (al), Since 7(2) obevs (a) to (e) of lemma 9, 


Nirv } min Nirwr*) 
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which (since the game is symmetric) implies that 
r(z) is optimal. Thus any optimal rf(z) must satisfy 


KNirw®)— min K(r,r*) 


r 


and thus, by lemma 1 and (c) of lemma 9, must 
satisfy rf{(1)—1. Hence the assumption at the end 
of section 4 is correct, and so, by lemma 8, r(2) is 
the on/y optimal strategy. 


7. Discussion of a Limiting Case 


The limiting case 6= 1 will now be briefly discussed. 
(Reeall that 6>1 previously.) It is easily verified 
that the results of the previous sections remain 
valid (many of the proofs can be considerably simpli- 
fied), so that the game with b=1 has a unique opti- 
mal strategy which is obtained from that for 6>1 
(i.e., from (a) to (d) of lemma 8) by setting b= 1. 
Since g=0 for b=1, the unique optimal strategy 
is given by 


2 r(z)—al(a+]l on lO.(a—-1)/a). 
(b ry(z)=1/(a+1 on \Oa-Lliva 
c r,(z)=1 on (a-1)/a,t}. 


Upon changing 
laver Ry plaver S. ro. ra. ra. 
phar) pia) ay V ‘ 


to plaver 1, plaver 2, ps, po, py, a/b 
respectively, one finds that this strategy is precisely 
the unique optimal strategy of the continuous von 
Neumann poker game discussed in sections 19.4 
to 19.10 of [3|. To explain this, note that (2.2) and 
(2.3), with b= 1, vield 


(7.(2)—G,(2) R, 2| [r.(t) +’r3(t) |\dt <0 


for all strategies 7(2), which on the one hand shows 
that the first course of action is dominated, and on 
the other hand vields a quick proof" that r)(2)-0 
on {0,1} for all optimal r(z). Thus, without loss of 
optimal strategies, the game can be “reduced” by 
“striking out” the first course of action for both 
players. This reduced game is easily seen to be 
identical with the von Neumann game, except for 
some changes of language and a multiplicative factor 
b in the expected payoff function. Thus the von 
Neumann game appears as a reduced form of the 
limiting case b= 1 of the game solved in this paper. 


tifrn(2 0 for some 2’, then 7? >0; the formula and the domination eriterion 
vield ri(z 0 for all z, a contradiction 
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Moebius Function on the Lattice of Dense Subgraphs 


R. E. Nettleton and M. S. Green 


(August 3, 1959 


The 


Moebius function f, on the lattice of k-dense subgraphs of a connected graph, 


defined in a previous paper, is calculated for graphs G containing isthmuses and articulators. 
ff, evaluated for the null graph @ is shown to vanish if G contains an isthmus, while for any 


integer g there exist graphs containing articulators for which f;(@) 


=q. The “lattice of path 


sets’ joining a pair of points and the lattice of graphs ‘associated with G and a subgraph G’”’ 
are defined and the Moebius functions on these lattices are shown in certain cases to be 


related to fy 


“k-dense subgraph,” “‘isthmus,”’ and 
1 


The concepts 
“articulator” were defined in a previous paper [1 
where it was shown that the ‘-dense subgraphs of 





‘a given connected graph G, together with the null 


graph ¢, form a lattice. The Moebius function /, on 
this lattice will be defined in the present paper, and 
f(¢) and, in particular, f,\(¢@) will be evaluated for 
various types of graphs. It is found that f\(@)=0 
if @ contains an isthmus, and | f\(¢)/<1 unless @ 
contains an articulator. However, for every  in- 
teger g, there exists a graph containing an articulator 
for which f\(@)=q. A second lattice is formed by 
the sets of paths joining a given pair of points and 
containing among them all the points of G, and when 
( contains an articulator in which there are just 
iwo points, f\(@) can be expressed in terms of the 
Moebius function defined on the second lattice. 
The graphs “associated with a graph G@ and a sub- 
graph G’” are defined and shown to form a lattice. 
A relation is found between the Moebius function 
on this lattice and the functions ti() defined for 
the individual lattice elements. 


1. Definitions 


A linear, undirected graph G is a set of elements, 
called points, together with a set of ordered pairs 
of these elements which define a symmetric non- 
reflexive, binary relation. For any given graph G, 
the number of points will be assumed finite and unless 
otherwise specified this number will be denoted by 
the symbol xn. A subset of the set of points in G, 
together with all the ordered pairs which contain 
only points in the subset, is called a subgraph. The 
union, intersection, or difference of two subgraphs is 
the subgraph determined respectively by the union, 
intersection or difference of the sets of points in the 
two subgraphs. If @ is a subgraph of G@, we say 
( is contained in Gand write G’ CG. Two points 
appearing in the same ordered pair in G@ are said to 


Figures in brackets indicate the literature references at the end of this paper 
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For any given pair of points, p and gq, 
a sequence of distinct points ¢{p, +» Pr}, with 
the property that py=p, p,—q, and p, and p,,, are 
neighbors for 1<7<r, is called a path joining p and gq, 
in which p,; is the tnitial point and p, the terminal 
point. A path containing r points will be called an 
r-cycle if its initial and terminal points are neighbors 
and no point of the path is a neighbor of more than 
two other points of the path. If the sequence of 
points of a path /? contains a proper subsequence 
which is also a path, the subsequence is called u 
subpath contained in ?. Two points joined by a 
path are said to be connected. A graph is connected 
if it has only one point, or more than one point and 
each pair of points connected. If every pair of points 
are neighbors, @ is said to be completely connected. 
A graph which not connected is disconnected. The 
null graph @ is disconnected. Unless otherwise 
specified, the symbol G will denote a connected 
graph 

If G’ is a subgraph and G-G’ is not connected, @’ 
is said to disconnect G. If G@’ disconnects G@ and 
(@-G’ contains, and is contained in, a set of connected 
graphs, the union of any pair of which is not con- 
nected, this set is called the partition of G-G’. If @’ 
contains precisely m points, the partition of G-G’ 
consists of at least k+1 connected graphs, and @’ 
contains no proper subgraph which disconnects G, 
then @’ is called an |m, k|-isthmus if it is completely 
connected and an [m, k]-articulator if it is not com- 
pletely connected. An [m, 1]-articulator or [m, 1]- 
isthmus will be called an articulator or isthmus 
respectively when the number of points is not 
relevant. 

For any subgraph G@’, G(@’) will denote the sub- 
graph determined by the set of all points in G-G”’ 
which have neighbors in G’. If @’ is a single point 
p, we shall denote this set by G@(p). A connected 
subgraph G@’ is said to be k-dense provided there are 
at most k-1 points in G-G@’ which are not points of 
G(G@’). A k-dense subgraph which contains no k- 
dense proper subgraph is [),-minima/l. A k-dense 


be neighbors. 








proper subgraph of G, which is contained in no other 
k-dense subgraph except @ itself, is called J),- 
maximal. The symbol [, will denote the union of 
all the ),-minimal subgraphs of G. 

Let S, denote the set of all the k-dense subgraphs 
of G which contain at least & points, together with 
the null graph. It has been —e in a previous 
paper [{1, theorem 2.2] that the set S, form a lattice 
under the relation of set inclusion, in whie h the lub. 
of two subgraphs is their union; and the g.l.b. of two 
subgraphs is their intersection, if it is in S;, and 
otherwise is either @ or a graph in the partition of 
the intersection, if the latter is not connected. The 
lattice formed by the graphs in S; will be called the 
lattice of k-dense subgraphs. For any finite lattice 
composed of elements partially ordered by the rela- 
tion > there is a greatest element / such that J/>,r 
for every element « [2, ch. IL]. The Moebius func- 
tion on the lattice is a relation which associates with 
every lattice element x a unique integral number 
M(x), defined by 


M(1)=1 
M(x) >My) 


y> 


For any lattice element 2, M(x) given by this 
definition is equal to the number given for the same 
element by the usual definition of the Moebius func- 
tion [2, ch. I] on the dual lattice obtained by replac- 
ing > by <. On the lattice of subgraphs in S, for 
which @ is the greatest element and = is understood 
to mean >, the Moebius function will be denoted 
by f, and the number which it associates with a 
given subgraph G@’ by f, (@’). A second function, 
defined on the entire set of non-null subgraphs of @ 
is the size o. For any subgraph G@’ #@, o(@’) is the 
number of points in @’. 


2. Completely Connected Graphs and 
Graphs Containing Isthmuses 


The Moebius function on the lattice of /-dense 
subgraphs will be shown to depend on the size of @ 
and on the way in which the points are connected. 
First, we shall prove some preliminary combinatorial! 
theorems. 

Lemma 2.1. If G,is a subgraph containing a proper 
l-dense subgraph G, and S is the set of all subgraphs 
of G, which are 1-dense and contain G@,, then? 


> (—)’ 


Ges 


, 
- 0). 


It has been proved [3, lemma 2.1] that a subgraph 
which contains a l-dense subgraph is 1-dense. 
Therefore a subgraph in S is formed by taking G, 
alone or the union of G, and the subgraph determined 


by any set of points from G,-G,. If ¢(G,) —¢(@,) =m, 


2? The notation X4.2 will denote the sum over all subgraphs .1 belonging to 
the set B of subgraphs, while SecA.” will denote the sum over all subgraphs 
1 which belong to B and which contain the subgraph C, and Zeros 4.4 the sum 
over all A in B which are contained in ¢ 
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the number of subgraphs in G,-G, 


' . (m 
points Is ( ) and so 


I 


X ( ¢ \ Mm k ‘ 
ZA) (—) ( l. )i yk — (—)o'@2)(] _] mig 
Ges 0 ’ 


A sequence of subgraphs |H,, . . /1,}, for any 
positive integer qg, will be « -alled an 1 H-sequenee if, 
and only if, it has the following properties: (1) H/, js 
connected; (2) for every 771, Hl, -G(H 

DVuiG(H) —H,. We W=—G’, for a given con- 
nected subgraph G’, the //-sequence is said to be 
hased on _ 


LemM, 2.2. lf HT], ; HT, is ait I1-sequence, 
then: 
For any r>1, >3. 47, is connected, 
lf q7 | and » & H1 is I-dense, then eithe 
there exists a non-null subgraph H,., such that 
H,. : /T,, HT, is an EH-se quence, or >»53 iH 


is I-dense. 

To prove part i suppose =p, is a point in H, for 
k<r. Then p, has a a Pe mm FT), and 
given a a p, in H, for +<k—1, p, has a neighbor 
pin IT, Thus there exis a path a= {eo 
p,} joining _ in //, to a point p; in 77; and whic ‘h 
contains no other point in //;. Similarly, if p’ is 
any point in >0;.,/7; distinct from p,, p’ is in some 
/1, for ar, and the ‘re exists a path /?” joining p’ to 
a point p,’ in /7; and containing no other point in 
H,. Since /[7/, is connec ‘ted, there is a path P, in H, 
joining p, and p,’. If these three paths are all dis- 


joint, i.e., no two have any point in common except 
for the distinct points p, and p,’, they form, take ‘n 
together, a path | p,, . + er ae ip 


joming p, and p’. If these paths are not disjoint 
except for distinet p,; and p,’ there exists a point p, 
in P77, which is also a point of P’. If 6 is the 


greatest integer for which this is the case, then the 


segments {p,,..., Po} and {p’,..., po}, of F 
and p’ respectively, when taken together, form a 
path {py, . .., p’) joining p, and p’. 


To prove part 2, suppose that there is no non-null 
H,., such that | 77, H1,, H1,.,) isan H-sequence. 


Then every point in G(//,) has a neighbor in one of 


the 77, for 1<i<g. Furthermore every point in H, 
has a neighbor in //, ,. From part 1, we know that 
p i/7, is connected and since it contains a neighbor | 


of every point in //7, and of every point which has a 
neighbor in >)%_,/7,, it must also be 1-dense. 


Lemma 2.3. Jf @’’ is a connected subgraph and @ 
is a I-dense subgraph containing G’’, then there is a 
unique H-sequence based on G’’ such that G’ is equal to 
the union of all the subgraphs in the H-sequence. 

Define 77,=@"’ and for 7 >1 let H,=@Q’'NGU1 
-S3'-3GU1,)—H,. Suppose that for some r, there 
are points of @ not in 7, At least one of 
these points, p say, must have a neighbor in some 
H, for some 7<r since otherwise @ would not be 
connected. If j<r, p must be in //,; for 7<r. This 
is contrary to hypothesis, and so p is in //,,,. Thus 
since has only a finite number of points, there 
must be a positive integer q such that a> JT, 


which have | 


a Se 


~~, A: ete 4 °F 


have } 


1)” () 


for any 
Pnee if. 
) A, is 
(HH, 
1 COon- 


UENCE, 
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h that 
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/], for 
1, and 
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WwW hie ‘h 
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is the 
n the 
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rhbor | 


has a 


nd (’ 
is a 
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Hy, 
there 
ne of 
some 

be 
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Thus 
there 


J, 


and it follows that G’ » a! since the //7, contain 
only points of G@’. Furthermore the sequence 
a _H,}, by its construction, is an //-sequence. 


Suppose that there is a second //-sequence 
Hi,. . ., Hi} based on @”’ such that @’- p> FAA 
H,=H/ for i<j and H,4H;. It then follows that 
H, must be unequal to G@’NG@(U1,_,)— D2i-i@ 
H, \— WH, =H. We have //i}c//,, and therefore 
there are points in /7,—/H; which are in graphs Me 
fork>j. Since these points have neighbors in /7)_, 
this cma is contrary to the de finition of an H-se- 
quence and thus H, IT;. Since /1, HT; G’’, it 
follows by induction that the two //-sequences are 
identical. 

An H-sequence based on a connected subgraph 
G’’ with the property that the union of all its graphs 
is 1-dense will be called a D-minimal //-sequence if it 
contains no proper subsequence which ts also an 
H- sequence based on @’’ with this property. By 2. 

} for a given connected subgraph Gr’ in G, to every 
1-dense subgraph G’ containing G’’ there corresponds 
a unique H]-sequence based on G@’’ such that @’ is 
the union of the graphs in the //-sequence. This 
[]-sequence contains a subsequence which is a D- 
minimal //-sequence based on G’’ and which is also 
unique. We can see this by observing that given 
two H-sequences {/7,’) and (/7,/’: which are sub- 
sequences of an J//-sequence (//,., if for some r, 
Hj=H =H, for all gsr and I),,A4Hy,, then 
either 17/., or 77,4), is aw ., and one of these must 
be equal to //7, for 7 >r+-1, which is impossible by 
definition of an + ile Thus (/7;, and (177 
must be identical if they are both J-minimal.  Ae- 
cordingly the /7/-sequences based on G’’ and corre- 
sponding to 1-dense subgraphs containing G’’ may 
be divided into families such that all the H-seque nees 
in each family contain a particular //-minimal //- 
sequence based on @’’, and then the //-sequence 
corresponding to a particular 1-dense subgraph con- 
taining @’’ will belong to one, and only one, family. 
Let S be the set of all l-dense subgraphs for which 
the corresponding //-sequences belong to a particular 
family; let (77,, . . ., HZ,) be the D-minimal //-se- 
quence contained in every //-sequence of the family; 
and let G,=Z4_/7,. Then every //-sequence in the 
family may be denoted by (//), .. ., Hy, Hy 
where /7,,, is the graph determined by any set of 
points, empty or nonempty, in G(/7,) —334-/G(1))- 
H1,=G,. There can be no //, for k>q+1 in any 
H-sequence of the family, since otherwise //, would 
have a neighbor in some //, for i<q, which is im- 
possible. By part 2 of 2.2, we know that G, is non- 
null if g>1. If g=1, G, is non-null if G’=—//, is a 
proper subgraph. Assume, therefore that G@’’ 4G, 
so that G, ~. The set S is the set of all 1-dense 
subgraphs in G,UG, which contain G@,. By 2.1, 
Por s(—)"°°? =0. and this must be true of every 
family of 7/-sequences based on G’’. Thus if S, 1 Is 
the set of all 1-dense subgraphs which contain @”’ 
 - tesa FO —() This establishes the theorem: 

THEOR -EM 2.4. If ’ wa proper connected sub- 
graph in G and Si is the set of all 1-dense subgraphs 
which contain G’ . then > a es, ( \a (Gl 0). 


—— 


The preliminary theorems and lemmas are now 
established which will make possible the determina- 
tion of the value of f, corresponding to any particular 
lattice point. Since the Moebius function is uniquely 
defined, any function must be equal to f, if it satisfies 
the equations which define f, rec ‘ursively on the 
lattice of k-dense subgraphs. Consider in particular 
the function g, whose domain is the set of subgraphs 
in S,, defined by 


gx.(G’) =(—)*te@) 


for 6#G’eS,NS,, 
gy G’)=0 


if G’ is in S,—@ but not 1-dense, and 


n(@)=(—)**!' SO (—-)@ 


o¥G'eSiqSe 


Since every subgraph which contains a 1-dense sub- 

graph is l-dense and thus k-dense, the set S of all 

subgraphs in S, which contain a given 1l-dense sub- 
'graph G@’’ in S, is the set of all subgraphs which 
contain @’’. Thus 


Sse") =(—)*4[ Sa(— ee —(— yr” | 


CG" <G'eSi— CG" Ges 


= )n+0(G"") n(G’’) : 


which follows by 2.1 if we set G;=G@ and G,=@"’. 
If @’’eS,—S,, the set of all subgraphs which properly 
contain @’’ and for which g, does not vanish is iden- 
tical with the set of S of all 1-dense subgraphs con- 
taining G’’. Since G@”’ is not 1l-dense, G’’#G, and 
since @’’ is k-dense, it is connected. Therefore, by 
2.4, 


>> = ge(G’)=—(—) *t* 2 (—)" "= 0=— 4 (C"’’). 


G"’ <G’eSi—G" G'eS 


Thus it is established that g,(@’) satisfies the re- 
currence relations which define f, ‘for all G’#@ in 
S,, and so g, (@’)=f,(@’) for G’A@. It then follows 
from its definition that g, (¢)=f,(@). Accordingly, 
| we have proved 

THeoreM 2.5. On any lattice of k-dense subgraphs, 

| Oh Se. 

Since for any non-null @’ in S,, f,(@’) depends 
only on o(@’) and on whether or not G’eS,, attention 
will be given to f,(@), and in particular f,(@), for 

| which the value depends on whether @ is completely 
connected, or contains an articulator or isthmus, 
and on other properties of G. 


Tureorem 2.6. Tf G@ is completely connected, 
| fi(o) =(—)". 

If G is completely connected, every point is a 
| neighbor of all other points and is therefore 1-dense, 
as, in fact, is every non-null subgraph. Thus, for 


, n 
every k>1, @ contains (") 1-dense subgraphs hav- 


43 








ing & points each, so that by 2.5, 


1) ¢ t= (—)", 


a result which follows from the binomial expansion 
of (1—1)". 


LEMMA 2.7. If S is a set of subgraphs of G, F, for 
j=] 9. n the set of all families of subgraphs 
in S such that each family contains precisely i sub- 
graphs, for any partic ular family C, BCC) is the union 
of all the subgraphs in C, ‘and S’ is the set of all sub- 
graphs of G which belong to a given set R and which 
contain one of the subgraphs in S, then 


f.@)=(—)"*" 9 ( 


k=1 


n 
Be (re ae 2 ee _ (- (A) 
es i=1 CeFy BIC)CG"ER 


Consider a particular subgraph @’eS’ and suppose 
that G@’ contains precisely g subgraphs from S which 
form a family F,. The subgraph G@’ will occur 
once in the triple sum in (A) for every subfamily 


. . . , ° ’ . ( 
of the family F,, and since F, contains ts) sub- 


families of precisely & graphs each, the total contri- 
bution to the triple sum in (A) of the summands 
corresponding to @’ is 


\o G’ 


(2) (r= 


Since this result holds for every subgraph @’ in 
S’, the right and left members of (A) are equal. 
THEOREM 2.8. If G contains at least k disjoint 


subgraphs, each of which disconnects G and at least 
one of which is an isthmus, f,(o)=0. 


If G,, G, are k disjoint subgraphs each of 
which disconnects G, then every 1-dense subgraph 
contains at least & points, one ey each G 
(¢=1, _k). Otherwise for some 7, (1 “<j <k). ) 
G—G, would be l-dense and thus connected Thus 


S; S, nS, and f.(o)=gn(e) =gi(@) by 2 Suppose 
one of the G,, G, say, is completely seanened 
that every subgraph of G, is a connected proper 
subgraph of G@ If S is the of non-null sub- 
graphs of G, and S’ the set of l-dense subgraphs 
containing a graph from the set S, we have S’=S), 
and thus 


set 


)o eo" 


Ji(o) =(—)"" 


If (’is any subfamily of S and B(C) the union of the 
graphs in B(C) is connected because G, is com- 
pletely connected, and thus by 2.4, 


> (—)o(o") —(. 
BIC) ecG”eS 


Since this is true of every family C, we find from 2.7 


that g,)(¢)=0. 
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Let us now suppose that S denotes the set of 
minimal subgraphs which belong to S, f) §, 
i.e., the subgraphs with this property which have 
no proper subgraphs belonging to S,; f S;, 
R is the set of 1-dense subgraphs. Since every 
l-dense subgraph having at least k points contains 
a minimal subgraph with this property, 2.5 implies 
if S’ is the set of all _——— in FP? which contain 


one of the subgraphs in S, then 
"7 / n+1 a) s(G’ 
Ik(d) =(—) (—) 
Ges’ 


If C’is any family of subgraphs in S and BCC) is the 
union of the graphs in then BCC) is connected. 
since the union of two l-dense subgraphs is 1-dense 
(2, lemma 2.1] and therefore connected. If @ js 
not equal to the union of all the subgraphs in §. 
b(C’) is a proper connected subgraph. Then by 2.4, 


> > { . }o(@") __._. 


— 
B(C)cG’eR 


Since this holds for every family ©, 2.7 implies 


yy ( \a(a" 


Cell 
(res 


QO; 


so that f/, vanishes. Thus we have proved 
THroreM 2.9. If Gis not equal to the union of all 
the minimal subgraphs with the property that they are 
I-dense and contain at leasi k points, then f.() (). 
From 2.9 is is seen that #,(@)=0 unless G is equal 
to the union of the J\-minimal subgraphs. If 
G=T, and that /,-minimal subgraphs are mutually 
disjoint, it has been proved elsewhere [{3, theorem 2.6 
that G is completely connected, that, by 2.6, 
hi(od) = ( ee" 
Thus we have proved 
THeorem 2.10. Jf the D\-minimal subgraphs ar 
mutually disjoint,  f\(o@) <1. 





so 


It has been proved (3, theorem 2.4] that if mn >1,@) 


contains at least two 7,-maximal subgraphs. Sup- 
pose n> and there is a /)\-minimal subgraph 6’ 
which is contained in every )),|-maximal sugraph, and 
7’ is D\-minimal, G’’#G’. Then there is a point 
p in @—G@", and @’+(G—p) is P,-maximal. 
This impossible, since this /);-maximal subgraph 
does not contain G’. Thus T,=4G", 
proper subgraph, that by 2.9, Ai(o) 0. This 
result is Hall’s theorem [4| for the spec inl case of - 


sO 


lattice of 1l-dense subgraphs. It may be stated 1 
the form: 
THEOREM i Ae On the lattice of 1-de WNE sub- 


graphs of G, if n>1 and @ is not the gdb. of any se 
of D,-marimal subgraphs, then f\(o)=0. 

It has been previously proved [3] that either @ is 
completely connected, in which ¢ f\(o),=1 by 
2.6, or @ contains a disconnecting subgraph which 
in turn must contain an articulator or an isthmus. 


case 


If G contains an isthmus, /,(¢)=0 by 2.7, and thus 
we have proved 
THEOREM 2.12. | f\(@)) <1 unless G contains an 


articulator. 


while } 
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3. Graphs Containing Articulators 


A simple example of a graph containing an ar- 
ticulator is an r-cycle thy r>3. Since each point is 
connected to only two other points, the two neighbors 
of any given point in the r-cyele constitute a dis- 
connecting subgraph G’, which not connected, 
each point of which is connected to both graphs 
in the partition of G-G’, and which is therefore a 
(2, 1|-articulator. 

ToeoremM 3.1. Jf G@ contains an n-cyele, 
flo) =—1 for kAn—1, and f, (¢)=n—1. 

if G contains an n-c ve ‘le, every point is one of a 
sequence (Pi, .- ~ , Py} such that each point is a 
neighbor only of those which immediately precede 
and immediately follow in the sequence, except for 
pi and p, which are neighbors. We have from 2.5, 


Is 


then 


hile) =(-)" rl,(—)* -1. For any 7, Gp, is 1-dense 
since P=! py. > Par Pr - + +> Pi-r} 18 @ path 
— ting Pix aed D, 1, Pi. being take ‘n to mean p, 
if 1=1, whie ‘h are neighbors of p,, and any two other 


-p, are joined by a subpath contained in 


nl in @ 
Pim py. is connected since 


P. Similarly for any 7, G 


there exists a path {pj, .. .. p containing all 
the other points. On - er rv hand, if p, and p, 
are not neighbors and a<b, there exists p, and pg 
such that a<ce<b and either d>b or d<a. It is 
easily shown by induction that any path joining 
p, and py must contain p, or p,, and thus G—p,—p» is 
not connected. Thus for n>3, S;—@ consists of G 
plus n subgraphs G—p, (¢=1, . n), and n sub- 
graphs G—p;—p;,, so that by 2.5) fi(@) 

(—)"*! [n(—)"-?+n(—)" !4+(—)"| 1. Also for 
n>2 we have f,_1(@) = ( ntlin(—)"-'+(—)"|=n-1. 


All graphs considered hitherto have been such 
that fi(@) <1. Consideration of graphs containing 


articulators, however, will show that there exist 
graphs for which f,(@) assumes arbitrarily large 


positive or negative values. First we must prove 

LemMaA 3.2. If Geontains an articulator G@’ such 
that each point in G—G’ is a neighbor of every point 
of GQ’ then 


>> (-) 


—_— 
G>G'esS 


fi@)=(—)"*"[1 


lf 77 is any ee containing a point p’ in @’ 
and a point pin G , then pand p’ are neighbors; 
every point in G’ is a neighbor of p; and every 
point in G—G" is a neighbor of p’. Thus the sub- 
graph p+ p’ is 1-dense as is /7 which contains it. 
If a((r’) 


mt ; 
m, there are ( ) subgraphs contained 


I 

. , ‘ } r 
(’ having k& points each and ( ) subgraphs con- 
having + points each which can be 


vl Mm 
j ) ) l-dense subgraphs, 
' , 


each with & points in G’ and r points in G—G@’. If 
Sis the set of all graphs in S; which have points in 


tained in G@ 


; m 
combined to give 


4 
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fi@)=( 


both G’ and G—G’, we have 


. (’ s( — mM 
(—)oun >( x tS * )( | sl. 
Hes. rei \ PF 


l-dense subgraphs contained in 
every 
2.5, 


Since there are no 
G—G’, which would otherwise be connected, 
graph in S,—Sj{ is contained in G’. Thus, by 


yy 


5 oe] 
G'>G"S:-—¢@ 
( efi . ye «| 


If q is any positive integer, there 
qand f\(¢) = —4q 


yuri )? H 


a | 


Tle 


. 


G'> G'S: -—¢@ 


THEOREM 3.3. 
exist graphs Ga nd G’ for which fio) 
respectively. 

Suppose there exists a graph G@, for which f,(@) =r. 
Consider the graph G, which has the following 
properties: G, consists of an articulator @)_ iso- 
morphic with G@, plus s points, p,, ps, at least 
two of which are not neighbors, and each of which 
is a neighbor of every point of G,. By 3.2, we have 
for G, the result that 

wie | 


Since any subgraph which is 1-dense in @), contains 


> 
2 iq; +8+ 
file) =(—)" | tt orate 


neighbors of p,, ..., p, and is therefore 1-dense 
in G,, it follows that 

f(b) =( \? G,)+1 — = \? Gq" 

fi Gi'> O%S,-6 
is the funetion f,(@) defined on the lattice of 1- 


dense subgraphs of G,. Furthermore, if two graphs 
are isomorphic, to every Il-dense subgraph of one 
there corresponds a unique l-dense subgraph of the 
other, and thus f,(@) for @, is equal to fi(d) for G@,, 


which we shall denote by f{(@). Thus 
fi() =(—) 7G FF + (—) *f 8 (@). 
If o(fG,) is even, f,(@)=(—)*{r—1}. In particular 


we can suppose G, contains an n-cycle, n even and 

2, so that r 1 and let « be odd, so that o(G,) 
isodd and f,(@) If s were even, then f,(@) = —2. 
However, we have also proved that if there exists 


= 


G, such that o(@,) is odd and f{(¢)=r, then there 
exists a graph G@, for which fi(¢@)=(—)'lr+1], 
containing o(G@,)+s points. If s is even, o(@) 1s 
odd and for Ge, fi(@=r+1. If is odd, fi(@) 
for Gis —[r+1}. By induction, it follows that 
for any integer g such that q 1 or |q|>2, there 
exists a graph G for which f,\(@)=q. Uf @ is com- 
pletely connected and n=2, fi(@)=1 by 2.6, and 


thus the theorem is proved for all positive integers q. 








THeoreM 3.4. If pand p’ are two distinct points of | 
G, and S is the set of all families of paths joining p and 
p’ such that each family contains n points among them, 
then the families of S, together with the null family, os, 
form a lattice under the relation of set inclusion. 

If a family F of paths contains all points of G, any 
family containing F is also in S, and thus the l.u.b. 
of two families in S is their union. If the intersec- 
tion of two families in S is not in S, it cannot con- 
tain any family of S. Thus the g.l.b. of two families 
is either their intersection or @s. For any pair of 
distinct points p and p’, the lattice formed by the 
families of paths which join p and p’, such that the 
paths of each family contain all points of G, will be 
called a lattice of path sets associated with G and join- 
ing pand p’. For any given lattice formed by a set 
of S of path sets joining a given pair of points, 7(/) 
will denote the number of distinct paths in each lat- 
tice element FA#¢s, two paths being called ‘“distinet”’ 
if they differ in at least one point, and ms will denote 
the greatest positive integer with the property that 
for some FY, mgs=r7(F). The symbol hs will denote 
the Moebius function on the lattice formed by path 
sets from S. 

THreoreM 3.5. If S is a set of families of paths, 
including the null family, which form a lattice of path 
sets associated w ith G. the n for each FeS, 


he(F’) =(—)™str 
for F#s, and if ms>1, 
he(és)=(—)™s#? 2 (—) 
os#leS 


The proof is very similar to that of 2.5 and will be 
left to the reader. 

If G contains a [2, k]-articulator consisting of points 
pand p’, every 1-dense subgraph must contain either 
por p’ since otherwise G— p— p’ would be 1-dense and 
therefore connected. Accordingly if RP is the set of 

1-dense subgraphs in @ and S’ the set of all subgraphs 


in FP which contain Pp, p’, or p+p’, we have R=S’, 
and then 2.5 and 2.7 imply that 

fi@=(—)**' So (oO = ex] >) (-—)e" 

Ces’ poG’eR 
‘ a(G ‘ 7(G"’ 
(J a A 

p’>G’'’eR p+p’oa"eR q 

The first two terms in the brackets vanish by 2.4 

since p and p’ are each connected proper sub- 


graphs. Every 1-dense subgraph containing both p 
and p’ must contain a path from the set S of all 
paths joining these two points, and so again em- 
ploying 2.7, we find that 


n 


fi(@)=(—)**' >) DD  (—ece’ 
i=1 CeF; B(C)CG"eR 
where F, is the set of all 7-ples of distinet 
paths from S, and (, B(C), are defined as in the 
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joining p and p’, 
£ / 


hypothesis of 2.7. For each family C, BCC) is cop. 
nected since the union of two connected subgraphs 
having a point in common is a 


ingly, by 2.4, >) (—)*‘””=0 unless B(C)=@ 
B(C)cG"eR ’ 
i.e., unless Cis an element of the lattice of path sets 


)7 Grr 


in which case >* ( 
B(C)CO"eR 


| (—)"*r©., Tf 8S” is the set of lattice elements, then 
Ai(¢) >) (—)'”, provided gs is not the 
8" FeS” 
only element in S’’, and f;(@)—0 otherwise. This 
result implies 
THeorem 3.6. Jf G@ contains a (2, k}-articulator 


and S is the set of elements Jorming the lattice of path 
sets associated with G and joining the points of the 
articulator, then 1\(¢)=(—)"shs(os) tf S contains? ‘a 
nonempty set of paths, and f.(6) 0 otherwise. 


Theorem 3.6 can be used to show that if @ con 
tains a [2,4]-articulator consisting of points p and p’ 
and f,;(@) #0, then f,(¢) ean be written as a product 
of factors, one factor for each connected graph in 
the partition of G-G’. Any path connecting p and 
p’ cannot contain points from more than one graph 
in the partition of G-G’, and so any set of paths 
forming one of the set S of elements of the lattice “of 
path sets joining p and p’ must contain a subset con- 
tained in G;+p-+ p’ and containing all the points of 
G,, for each connected graph G; (@=1, . , w> 
k+-1) in the partition of G-G’. Thus if S; ((=1, 

, w) denotes the elements of the lattice of path 


sets associated with G,+ p+ p’ and joining p and p’, 
we have by 3.5, 

u“” w 

m l ‘ * hdd l 

hs(@s)=(—)" 11 p > yrit ) (—)™st+! 

i=l \.os #G@G' eS ] 

{ s he (Dx 

Since ms=>\ms., we have established 


THroreM 3.7. Lf Geontainsa [2, k\-articulator G’, 
G,; for i=1,...,k+1 are the ead graphs in 
the partition of G-G’, and for each G;, 8S; 
elements which form the lattice of path sets associated 
with G,OG’ and joining the points of the articulator, 


| then filo) X() im plic s 


k+] 
4A (0) ; Il h. (® 
i=] U 


4. Associated Graphs 


Given a graph G@ and a subgraph G’ which is not 
properly contained in a completely connected ny 
graph, the set of graphs associated with Gand @’ 
defined to be set of all graphs @’’ which have the 
following properties: (1) A one-one correspondence 
exists between the points of Gand those of @’’ such 


is the set ot 


Accord. 


- 


is CO}. 
graphs 


Lccord. | 


(’) (i, 
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3, then 
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of the 


Wns ta 
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ind p’ 
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that neighbors in G are mapped into neighbors in 
Gi’: (2) @”’ contains a maximal completely connected 
subgraph, which will be called A[@’’}, with the prop- 
erty that A[@’’| contains the points in a” which cor- 
respond to @ in G; (3) any pair of neighbors in @’”’ 
such that both are not in A|@’’], is mapped into a 
pair of neighbors in G. If G, and G, are any two 
graphs in the set associated with G and G’, we shall 
sav G,>G, whenever the points in @ corresponding 
to AJG,| are a subset of the points corresponding to 
K1G,|, and G, >G@, whenever G,>G, but G,#@). 

THeoreM 4.1. Given G and a subgraph G’, the set 
of graph s associated with Gand G’, if nonempty, form 
a lattice under the relation >. 

The proof will be left to the reader. We shall 
denote by h(@’’; G, G’) the Moebius function evalu- 
ated for any G@’’ in the lattice of graphs associated 
with @ and G’, while f\(@; @’’) will denote the 
Moebius function evaluated for @ on the lattice of 
subgraphs I-dense in G’’. The function Ah, whose 
domain is the set of graphs associated with @ and 
@’ is defined by h,(@’; @, @’) =(— rv *'@"D +", where 
G’’ is any graph of the set. 

TueoreM 4.2. On the lattice formed by the set of 
graphs associated with G and Gi’. h, h. 

If G,, is the greatest element of the lattice A[G,,| 

(1, and o(KIG,,|) ”, so that h(G,: G, fr’) # If 
Qr< G,, and w o( A{G,,|) o( K{G’’)) there are (?) 
graphs /7 which satisfy the condition that @’’: H 
<i, and o( A|/7|)—eT1[G"’|\)=k. Therefore 


> > h,(H:;G, G’)=(—)e hte tet 
H>G@” 
= Ww ” ‘ 
x 21, ) &— (—)e X10 +9—5, (Q’’; G,G@’). 
k=1 - 


Since A, satisfies the recurrence relations which 
uniquely determine h, the theorem follows. 

Let R{G’’| be the set of subgraphs 1-dense in @’’ 
and Q/G’’| the set of all subgraphs in R[G@’’| which 
contain at least one point from A[G@’’}. 

LemMa 4.3. Jf G’’ is any graph in the set associated 


with Gand G’ such that K|G"| 4G", then 


AGE %)=(—)O""" 3 (—)yr™. 
HeR{(G'—Q{e”) 
by 2.5, 
f, (o;@’’)=(—)o Pt! 2, (—)@ 
17eQ(G") 
+ : (—)¢ | 
MeR(GT-Q@")} 

By 2.7, the first sum in the bracket can be expressed 


as a sum over families of subsets of A[@’’|. If Cis 
such a family and B(C) the union of the graphs in 
the family, B(C) is a connected proper subgraph 
because A[@’’| is completely connected and #@”’, 
and thus 

» S ( 


a 


Bic ‘> Hek(@") 


\o H). 0 


by 2.4. Since this is true for every family C 
subgraphs of A|@’’|, we have 
S* (— 


He@a’”} 


ait) 0 


TueoreMm 4.4. If @’’ is a graph belonging to the 


| set associated with Gand G’, then 


K|G’’| is 1-dense in G’’, 


>> AUT; G,G)f,(o; H)=0. 


H>G" 


1) 1G” 


2) If G’—KlG@"’| is not 1-dense in G”’, 


>> h(H; G, G@’)f,(¢; H) =(—)"@" 
1>Ga" 


I 
By 4.2 and 4.3, we have for G’’ 4#AK[G@’’| 


>> > 
—! 

H>G” = H'eR(H)—Q[H] 

( yo (KH) )+o(H’) 4 


>>) AUT: G, GA’) Sf (e, ID 


H>G" 
)” +1 


(B) 


If //’ is any subgraph of G’’ which is one of the 


_ set S associated with G and @’, 7H’ corresponds to a 


| [3] R. EB. 
of | 
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unique subgraph 77; CG. In turn 77, corresponds 
to a unique subgraph /7’’ in any other graph of 
S, and these two correspondences determine a 
correspondence between i and H’’. If H>@”, 
every subgraph of //—A{//| corresponds to a sub- 
graph in G’’—K|@’’|. If @’’ — K[@"’’| is not 1-dense, 
it contains no 1-dense subgraph, and so_ there 
can be no corresponding 1-dense subgraph in 
H1—K\H) and Ri7\|—Q|H) is empty for > GQ’. 
Therefore, if @’ AAK[G’’|, —h(G”: G, G’) file: @”’) 
can be added to both members of (B), and then the 
double sum vanishes, proving part 2 for this case. 
lf G’’= K{@"’|, filo; @’’) = (— )" and part 2 is obvious. 
Suppose G’’—K{[@’’] is 1-dense, and H’eR[@’’ 
Q(G’|. If w=o(G”"’ —H’—K|@"’ |) andk=1,. . ., 
w, there are (7) graphs /7>G@"’ with the property 
that A|//|=A[G"’|+k and H—K{|H) contains a 
subgraph which corresponds to 7’. The contribu- 
tion to the double sum in (B) of all the graphs which 
correspond to a particular /7’eR(@’’}|—Q[@’| is 


» 3 he 
f= \k 


) yo (K(G") +k+o(H") __(__ yo (K(@"]) +1 +0(H" 
1 


for 
Thus 


contains a similar contribution 
Q1G’’| except G’’ —K[@’]. 


(B) 
MW eR(G"| 


every 


p> hh; GQ’) f (6; DI) =hiG'"';: G,@ fo; @"") 
GG" 
(—)o(K1@""})+1+4+0(G"" 


K(G’"}) | ( etl. 


Since o(@’’)=0(G@)=n, part 1 follows. 
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The Minimum of a Certain Linear Form’ 


Karl Goldberg 


(Mareh 23, 1959) 


The positive minimum of the integral linear form L(a, . . 
Ti, . . “* 


found subjeet to the conditions a; >0 and L( 


Let a, <a. < <a, be n>3 positive integers. 
We seek the positive minimum JV of the linear 
form 

L (2,22, gl, ) =A, 


T tS T T Cyl n 


over all non-negative integers Ba, Bis: tg Be such that 


Lay .fe, L,) > 2a, a, (1) 
for all ¢=1, 2, 

Let |a;,d,| « 
al; and (>. 

For each 7=3, 4, ..., m, define r, in the following 
wav: If either a, or a, divides a,, or if a; =a, for some 
jti, set r,=0. Otherwise, let 7; be the minimum 
of the least non-negative residues modulo a, of 


ver, - 
lenote the least common multiple of 


a,—a;, 2a,.—a;,.. ., |(a;—1)/az]a2—a,;. 


We shall prove 


THeorem: M is the minimum of 2{a,,a2\, 2a3+ 


ee 2a,+Ta- 


3, 


2s T 
As a consequence we have the inequality 


2a;-+-a,—1>M>2ay.. 


Also, if L(y... ., t,) =, then at most three of the 
“ are positive. At least two must be positive. 
exactly two are positive, then either 2,;=[a@;,@2\/a, 
and a=|d),d2|/d2, or 4,=a,/a, and £,—1, or 22=4;/dy 
and 4,=1, or #;=42;,—1 for some 7 >i>3. If three 
of the 4, are positive, then both 2, and x are positive; 
the other positive zs, equals 1 and we have 2, 
—~{(dt.—a,)/a,) for that 7 Under any conditions M 
is achieved only with 2,<1 for all 7>3. We shall 
prove all this. 

M. Newman ? refers to our theorem in the case 
3. Weshall treat this case first. 

We have a,<a.<a;, and we want to find the 


Mt 





! The preparation of this paper was supported in part by the Office of Naval 
Research. 

2M. Newman, Construction and application of a class of modular functions, 
Il, Proc. London Math. Soc. 9, 373 (1959). 


If | 
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is 


“* I») 
r,) >2a,7; fori 


ay,t) 
1,2 


Onl n 
» *. 


positive minimum .V/ of the linear form 
Lay 2,23) = yt, + Ayla + Agry 


over all non-negative integers 4), 22, 3 satisfying 


1. (ay ,dt2,t3) > 2a id’ 4 1 ie ea (2) 
Let x;=—|(a,—a,)/a,|. Beeause a,<a;, 2; is 
non-negative. It satisfies 
a,—1>a,.—a,+a,2;>0. (3) 
Because a, <dy, this implies 
3 > As (2 a,)— 1>a,z}. (4) 
Now consider L(xj,1,1)=a,aj+a2+a,. From (3) 
we have 
2a,-+4 a,—1 > L(x},1,1) > 2as. (5) 


We know that 2a,>2a,, so that L(z},1,1) > 2a, > 2a». 
Finally, (4) yields L(z},1,1) >2a, >2a,2;. 

This proves that 2,;=.2), 42.=2;=1 satisfies (2). It 
follows that the left-hand inequality in (5) holds 
for M: 


9 


— 


2a,;,+a,—1>M. (6) 
From this point we assume that 2, 2), 23 satisfy 
(2) and 


M. 


L(y 22,23) 


Since L(2,,22,03) >2asr,, we have from (6) and 


a,>a, that 7,=0 or r,=1. 

If x;=0, than (2) implies a,2;—a.%. Under this 
condition the minimum value of L(2,22,2%3) is 2[a),a@»), 
occurring for x2, =[d),@2]/a, and 2.=|a),d2]/d>. 

From now on z,=1. From M=a,x,+a.7,+ a, and 
(2), we have M>2a,. From (6) we have 


2a,+a,—1 > 2a3+ 4,2; + (dgr%.—Aas), 


from which it follows that z,;>0 implies a;—1 > a@y2». 
If 2, =0, then (2) implies a.27. = d73, so that M = 2{a»,a,). 








But 2{a,,a,] >2a;+a,—1 unless ay divides a3. Thus 
0 is possible only if a, divides as, in which case 
d3/da, and M=2a,. Similarly z,=0 is possible only 
if a, divides ay, in which case 2,;—a,/a, and Mf=2as. 
Since as divisible by either a, or dy leads to M= 2a, 
which is the best possible result with 2;—1, we may 
now assume that neither a, nor a, divides a, and that 
Lys >0. 

With 2, >0 we must have a4,—1>dyr. Fix 2. We 
shall find that permissible value of 2, which mini- 
mizes L(x ,a'2,23) ay; +r, +a. Clearly this is the 
least positive value of 2, satisfying (2). We have 


ry 


I's 


1. (ay 22,43) = Aya, + 1 + (ag — 1 — dary) + 2d ory >2a.2, 


for any value of «,. The other inequalities require 


a3 


Since 2dgr,> 2a, >a,, there are values of x, satisfying 
these inequalities. The least such x, is the least in- 
teger greater than or equal to (@3—dy22)/a,. This last 
quantity is positive, so this value of x, is positive. 
It can be written 


” M21 2— Az 
I) om » 
ay 


Let rs(22) be the least non-negative residue modulo 
a, Of dot.—a,. Then rg(ar.)=dor.—a,+a,r;'. It fol- 
lows that E 
L(zy’ .Xe,1)=2a3+13(2e). 
We want the least of these values for #: lving between 
1 and [(a;—1)/a,). Under our assumptions on the 
divisibility of a, this is just 2a;-+-7; with ry as defined 
in the theorem. This proves the theorem for n=3. 
Now assume 7 >3. We have a, <a.< s., 


and we want to find the positive minimum VV of the 


50 


linear form L(s,.25, r,) over all non-negatiye 
Integers J),f2, . . ., £, Satisfying (1). 
If x),42,%3 satisfy then » @.<<e 


(2), ; 
satisfy (1). Therefore our new M satisfies (6). Ley? 


Ly Lae 


L (2,2, ‘ 09 Sa) M. 
Then 
2a34+-@,— 1 Sa, h, + Gols +- A3%34 .. +a,7 
>a, (2, +22) +d3(43+ . +a,), 
It follows that 
F ry +F,<2 
and that s-+-44- .. +4,=2 requires s,+2,=0 


On the other hand, 


3a; >2a,+a,—1> L(a,,4 r,) > 2a, 
2 54, n, 
implies 
2 * %. 2 4 n 
Assume 2#3+-3,4+ . +e, =2. Then z,;=2;=] 
for some 1j>3 and all other s,—0. Then (1 


implies a;—a,; and \,=2a,. Again this is the best 


possible result with 2 l. 


If wtat+...+2,=0, then (1) implies Aydt 
dot). As before, this implies = 2{a,,a,). 
Ife,ta+... +2, 1, then x+,—1 for some 1>3 


and all other v,—0 for k>3. 
reverts to the case n=3 with a; replacing as. Our 
previous arguments complete the proof of the 
theorem, and the statements in the subsequent 
paragraph. 
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Space of k-Commutative Matrices’ 
Marvin Marcus? and N. A. Khan 


(August 14, 1959) 


Let [A, X]=AX 


mute” 
space of the linear transformation T(Y) 


XA and [A, X]}, 


1. Introduction 


Let A be a fixed N-square complex matrix and let 





[A, X]=AX—XA, [A, X)x=[A, [A, X]e-,]. It is 
easily checked that 
k k 
A,X);=>>5(—1)' AP-*XA*. (1.1 
1X h= 29 *() ) 
It is clear that the set of X such that [A, X],=0 is 


a linear subspace of the space My of all N-square 
complex matrices. This subspace is denoted by 
(i,(A). In theorem 1 is determined the dimension 
of (A) in terms of the degrees of the elementary 
divisors of A when there is exactly one elementary 
divisor for each eigenvalue. Let /, denote the set 
of matrices in My with precisely g distinct eigen- 


values. In theorem 2 it is shown that in case 
k>2(N—q) +1, then 
min dim (,(.A) = R(Q+ 1)?+ (¢q— R)& 
AcE gM» 
where N=Qq+R, O<R<q. The maximum is also 
found. 
2. Results 


Let 7 denote the linear transformation on My de- 
fined by T(.X)=[A, X] and we note that 7*CX) 
(A, X],. With respect to the basis /;, in My, or- 
dered lexicographically, we check that 7 has the 
matrix representation Jy@A—A’@ly. The nota- 
tion is the following: /,, is the N-square matrix with 
1 in position 7, 7, 0 elsewhere; /y@A denotes the 
Kronecker product of the N-square identity matrix 
ly with A. It is clear that one may assume A is in 
Jordan canonical form 


(2.1) 


s=1 


A 


where S* indicates direct sum and the J, are the 
Jordan blocks corresponding to the distinct eigen- 


!'This work was supported in part by U.S. National Science Foundation 
Grant NSFG-5416. z f 
2 Present address: The University of British Columbia, Vancouver, Canada, 


3 Present address: Muslim University, Aligarh, India. 


(A, [A, X] 4-1]. 


with a fixed matrix A are investigated. 





| 
51 


Those matrices X which ‘k-com- 
In particular, the dimension of the null 


[A, X], when A is nonderogatory is determined. 


values X,, s=1, ,qg, of A. If 7.(X)=0 and 
X is partitioned conformally with the partitioning 
(2.1) of A, Roth * shows that 


X=>°° X, 


s=l 


where X, is the same sizeasJ,, s=1, ,q. If, 
further, each J, is decomposed into a direct sum of 
companion matrices of the elementary divisors cor- 
responding to A, one may also effect a conformal 
partitioning of the corresponding X,. It is also clear 
that one may take A,=0 in examining the structure 
of X, since [J,, X,], remains invariant upon transla- 
tion of J,. We are thus reduced to considering the 
following situation in determining the dimension of 
(A) ° Let 

U lA 


t 


r;> 2 >ryp 


i=l 
where U,, is an r-square auxiliary unit matrix, an 
unbroken line of 1’s along the first super-diagonal 
0’s elsewhere, and suppose [U’, Y],=0. Partition Y 


conformally with U, Y=(¥,,) i,7=1,. . . ,l where 
Y,, is r;Xr;, and by (1.1) 
: 2) _ 
>> (—1y"({)e VU =0 i,j=1,...,1 
s=0 " 
2.2) 


is equivalent to [U’, Y|,=0. The problem then is to 
determine the number of arbitrary parameters in 
each Y4;. 

Equation (2.2) for a fixed 7, 7 represents a linear 
transformation mapping )’,,; into 0, and with respect 
to a suitably chosen basis this transformation has the 
matrix representation 74, where 

T7;= (1,,8U,,—U,,@1,,). (2.3) 
To simplify the notation, put r;=n, r,=m where it 
can be assumed without loss of generality that n>m. 
The similarity invariants of 7;, as computed by 


‘W. E. Roth, On k-commutative matrices, Trans, Am. Math. Soc. 39, 483 
(1936) 








Roth ° are f,(z) Cele h, fie-nsel(t)* 

gitte-| »y=—1, ,m and A=n—m. Hence 77; 
is similar to the direct sum of thee ompanion matrices | 
of these nontrivial similarity invariants, | 


m™m 
Ty >5° Clrot??-'), (2.4) 
p=1 


The sizes of these companion matrices arranged in 


decreasing order are A+2m—1, A+2m—3, 
A+3, A+1. 
Now, if k>A+2p—1, then (C(az 4*°?-"))*¥=0. 
If kc A+2p—1, then 
p({ C(aat2”-!) }*) —A+2n—1—k, | 
where p denotes rank. | 
Let » denote nullity. 
Lemmal: (a) 9(T%,)—km if 1<k<A, | 
‘rye k -A . 
(b) n(Ti;)=km (= = ) +C | 
if A<k<m+n—1, | 
(c) (Tt) =mn if k>m+n—1, | 


where C is 0 or 1/4 according as (k—A) ts even or odd. 


PRooF: 
(a) 1<k<. 
Then 
p(T%;) pa +-2p—1—k)=mn-mk, 
and n( Tt;) =mk. 
(b) A<k<n+m—1. 


Assume k— Ais odd and observe that the size of the | 
(m—(k—A-+-1)/2)th companion matrix in (2.4) 


A+-2m ( 2( nn 


Is 


A 


> )-1) k+2, 


and the size of the next companion matrix is k. 


Hence, 


p (T%,)=(A+2m—1—k) + (A+2m—3—k)-4 


rank of the kth 


)th companion matrix 


term in this sum is the 
-A+1)/2 


The last 
power of the (m—(k 


§'W. E. Roth, On direct product matrices, Bull. Am. Math. Soc. 40, 461 (1934). 
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| in (2.4). Thus, 
; F F 0! 
_ 2m—k-+-A—1\f/2m—k+A4+ 1 il 
p(Ti,) =( > \(: 5 ) ,! 
(AA— A) l 
(m—2=0y 1 
“ ; 
In case k—A is even, it is observed that the size of , 
the (m—(k—A)/2)th companion matrix in (2.4) js , 
k+-1. Also the size of the next companion matriy 
isk—1. Hence, 
p(T%;) (A+2m—1—k)+ (At+2m—3—k)- W 
\ 
r(k+1-b i 
(k—A)\? 
( m— > ) ° 
Hence, in either case 
n ( Ti) mn— p(T) ( 
(k—A)\? , ,, j 
mn—( R———— ) +0 ) @ 
r t 
k—A\’ , ,, 
mk—( > ) +( 
where C is 0 or 1/4 depending on whether k—A is 
even or odd. 
(c) k>m+n—1. 
Then 
C(x 2m De 0 
and 
n(T},)=mn. 
THeoreM 1. Assume A is N-square with distinc! 
ele nvalues as a ks Ay and let (x A,)* he the elemen- 
re divisors of A, j=1,..., q, 426@>...26 
Partition the integers 1,..., q so that 
€) €4, "Eo, 
€, Ces. | Ey, =—& ¢ 
Then 


(2) dim (,(A) kN—q(“—c) of k<2e,—1, 


aha e -_ ll 

(77) dim A,(A)=k Doe + SD 6-4 ~~ 
j=1 Qe-i tl , 

fe 1<}- 2€q, ' l 


Py a ° y 
(271i) dim (,(A)=>) é 1, where Cw 


j=1 


af k> 2; 


0 or % according as k is even or odd. 


' size of 
(2.4) js 
Matrix 


k—A is 


listinel 


rlemen- 
2 


ne fy is 


Since there is only one elementary divisor 


PROOF: { ) 
of A for each eigenvalue of A it may be assumed, as 
in (2.1), that 
a 
A . Jas 
s=1 
where J.= rd +U., UO, the e-square auxiliary 
ynit matrix, s=1, ..., g. Then, if XY is contained 
in (1,(1), 
. Se r 
X=) °° X, 
s=l 
where X, is ¢-square, s=1,..., g. By lemma 1 


we know that the number of arbitrary parameters 
in VY, is (by putting m=n=e,) 


i ‘ = 
h) Ne oe +-( if hk. 2e, ‘. 
4 
(¢) € if k > 2e, l. 
Consider (72) first: k-< 2e,—1. Then k- 2¢,,—1 
j=1,..., /, and hence XX, has n, arbitrary param- 
eters in it. There are gs—qe—; values of s such 
that ns=",. (qo=O for convenience). Hence, 
. a l 
dim (y.(A)=35 ns=D5 (qe—o-1)%q, 
s=1 a=] 
i Je 
yy Pen. ’ 
2a (Yo Yo )( ke até ) 
: } 
* _ aces 
2 (de Yo )€ a(", ( ) 
a | 
hk > €; a( —/( ) 
. 
kN a( a ): 
Next assume that (11) 2€,, 1<k- 2€,, ak. 
In this case 3 t / and hk e C 
MS CASC Ny = €],, 9 »-- yl aNd Ny €ig— gt ‘ 
o=1....,¢—1l. Hence, 
=! 
dim (,(A) 24 (de Je ny +> (Yo— Yo de. 
o— o=l 
1 a je 
mtSia+  a-0.(£-c) 
s—1 3 Qy—, +1 4 
If (i) k>2Qe,—1, then k>2e,—1 for j=1, +4 
and 
. , q 9 
dim Q,(A)=S3é. 
s=1 


In case there is more than one elementary divisor 
corresponding to a particular eigenvalue there does 


| not seem to be any simple formula for dim (,(A) 
in terms of the degrees of the elementary divisors of 


by repeated use of lemma 1, it is pos- 
(A) for any partic ular A. 


A. However, 
sible to compute dim (,( 
For example, if 








Ss t 6 
0 2 l 0 2 | 
A + T (2), 
0) 0) 2 1 0 2 
is ¢ 8 
then dim (,(A)=37 for k 
Next are determined the largest and smallest 
values that dim (,(A) may take on as A varies 


the set of matrices with precisely 4 distinct 


over E,, 
under the condition that k >2(J V—q)+1. 


eigenvalues, 


Lemma 2. Jf «> >e, are positive integers 
satisfying 
2 e,—N=qQ+R, 0<R<4q, 
then 
, > £. > 
R(Q+-1)?+ (q—-R)G< e? < (N—q+1)?+(q—1) 
j=l 
The lower bound is achieved for 
€) see €r > f €r+1 se €e Q, 
and the upper bound is achieved for 
e,=N—q+1 and e=... =e =1. 
Proor. The lower inequality is proved by induc- 
tion on R. In case R=O0, then: N/q=Q, and if 
€;,. . ., €, are regarded as continuous variables, then 
d . . . * 
3 @hasaminimum fore;=Q, j=1,...,9qg. Now 
j=1 
suppose the result is true for all remainders obtained 
by dividing N by q that are less than R. We first 
claim that there exists an integer iq such that 
€, >€:4, and €,>Q+1. Clearly the set of integers 7 
such that «,>Q+1 is nonempty otherwise (since 
R>~0), 
i. 
24 5S 90< N 
j= 
Let 7 be the largest integer 7 such that e,>Q+1; 
then if i were g, «>. >e,>Q+1 and 
Sey g(Q+y>N. 
j=l 
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Hence i<q and from the definition of i,  €:>e€:4;. 
Let My €;, JF and Me=e;—1. Then My > —— 
Hie Mt... tae=N—-1=Qq+R-1, and by 
induction 


(q—(R—1))Q@ 
+(q—R)Q 


R(Q | 1)? 
+ (Q?—(Q+1)’). 


Soyu3> (R—1)(Q-+1)24 
=I 


Now 


q 
ae , 
ej—2e,+1, 


j= j=1 
and thus, 


S> R(Q+1)?+ (q—R)Q?+2(e,—Q—1). 
j=1 


Since ¢,>Q-+-1, the proof is complete. 
bound is easily obtained. 
THEOREM 2. Jf k>2(N-q)+1, then 


The upper 


min dim (,(A)=R(Q+1)?+ (q—R)Q@ 


AcE My 
and 
max dim (,.(A)=(N—q+1)?+q—1, 
AcE Qh Mn 
where 


N= qQ +R, 0< R<q. 
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Proor. Let Ae“, and suppose «> .. . >e, ar 
such integers that ¢, is the sum of the degrees of glj 
elementary divisors of A corresponding to d, 
. ’ 


eek, 2s is 6 & 
Then 
q , 
2 t=, 
i 


— 


and hence, ¢ <N—gq+1 and 2e,—1<k. Thus kj 
at least 2u—1 where yu is the degree of any elemep. 
tary divisor of A. From lemma 1 one may cheek 
in this case that min dim (,(A) may be evaluated 
AcE gn Mwy 

by confining A to those matrices having precisely one 
elementary divisor for each eigenvalue. Hence 
(t—),)e, may be taken as the elementary divisors 
of A, gl, . ss q. By theorem 1 if AeZ,, 


. q ” 
dim (i, (A)=35 é, 
j=1 


and the results follow from lemma 2. 


Wasuineaton, D.C, (Paper 64B1-21) 


' 


a GQ Oe at a Gate om Eee 


> €, are 
'S of all 
to hs, 


us & jg 
‘lemen- 
* check 
nluated 


ely one 


Hence 
livisors 


31-21) 


JOURNAL OF RESEARCH of the National Bureau of Standa 






s—B. Mathematics and Mathematical Physics 


Vol. 64B, No. 1, January-March 1960 


5 


Selected Bibliography of Statistical Literature, 1930 to 
: 1957: I. Correlation and Regression Theory 


Lola S. Deming 


' 


(October 16, 1959) 


This is the first in a series of bibliographies dealing with various specific subjects in 


the field of statistics. 


References and titles of important contributions in correlation and 


regression theory have been taken from technical journals published throughout the world 


since 1930. 


Statistics per se is a relatively new science. Its 
| literature has grown tremendously during the last 
several decades. Not the least of the statistician’s 
difficulties is keeping abreast of this ever increasing 
abundance of writing even in his own particular 
area of the subject. Acquisition of more than a 
superficial acquaintance with statistical develop- 
ments is rendered enormously more difficult by the 
fact that publications in this highly applied science 
appear in a myriad of journals, many of which are 
not readily accessible. This bibliography is in- 
tended to be of service by listing the most important 
contributions in statistical theory and methodology 
as judged by two prominent reviewing journals. 
The NBS Statistical Engineering Laboratory main- 
tains a card file of abstracts taken from Zentralblatt 
fiir Mathematik for the years 1930 to 1939 and from 
‘Mathematical Reviews from 1940 onward. Each 
abstract has been assigned at least one subject 
category following the classification scheme of 
Mathematical Reviews. The present bibliography 
' lists publications relating to correlation and regres- 

sion theory. It is planned to follow this publication 
with similar ones dealing with other subject cate- 
gories of statistics. As later abstracts appear and 
are classified and processed, it is hoped that they too 
will be published as supplements to the present 
bibliography. A considerable time lag is unavoid- 


— 


’ able, however. 


, Teports, etc.) are italicized, and are followed by the | 


The following information is extracted directly 
from the abstracts: 

Author: The author’s surname is followed by 
initials only. In the case of complicated surnames, 
the first capitalized word given in the reviewing 
journal is used. When the paper bears the names of 
more than one author, the journal reference appears 
with each author’s name but the title of the paper 
is given with the senior author only. Multiple 
authorship is denoted by the symbol @ preceding 
the surname. 

Title: Given exactly as in the reviewing journal. 
Titles of separately bound publications (books, 


publisher. 
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Reference to literature: The references were tran- 
scribed to punched cards thereby necessitating severe 
and unconventional abbreviations. The name of 
the journal in italics and the number of the volume 
in bold face, are followed by initial page number. 

Date of publication: The next figure, in parentheses, 
shows the date when the article or book itself ap- 
peared. 

M (for Mathematical Reviews) and Z (for Zentral- 
blatt fiir Mathematik) are followed by the volume 
number and page number of the reviewing journal 
in which the abstract appears. 


It is a pleasure to record the help of several mem- 
bers of the Statistical Engineering Laboratory. 
Dr. Churchill Eisenhart, who originally conceived 
the plans for this project, has always been ready 
with assistance: and encouragement. Dr. Eugene 
Lukacs, now at the Catholic University of America, 
studied the alternatives for the subject classification 
and recommended the one used here. He also made 
innumerable selections from reviews not classified 
under the regular headings of probability and statis- 
tics. Dr. Norman Severo, now at the University of 
Buffalo, helped with the plans for mechanizing and 
transcribing the references in a form ouleallle for 
publication. 
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This is the second in a series of bibliographies dealing with various specific subjects 


in the field of statisties. 


References and titles of important contributions to the study of 


time series have been taken from a wide variety of technical journals published in the many 
languages and countries which have been actively engaged in statistical analysis. 


If this series of subject-classified bibliographies can 
assist the overburdened statistician to keep abreast 
of the great abundance of literature in his field, our 
aim has been accomplished. Although the coverage 
of technical journals is worldwide, each bibliography 
is still far from complete. The books and articles 
considered by two prominent reviewing journals to 
be of sufficient statistical importance to include in 
their regular abstracting services form the source 
material for these classified bibliographies. 

This particular subject classification on Time Series 
follows the earlier and closely related one on Correla- 
tion and Regression Theory.’ The titles and references 
are extracted from a card file made up of abstracts 
taken from Zentralblatt fiir Mathematik for the years 
1930 to 1939, and from Mathematical Reviews from 
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tained on a current basis in the NBS Statistical Engi- | 


neering Laboratory ; the bibliographies extracted from 
it are necessarily delayed due to loss of time in proc- 
essing, editing, and publishing. 

The following information comes directly from the 
abstracts: 

Author: The author’s surname is followed by ini- 
tials only. In the case of complicated surnames, the 
first capitalized word given in the reviewing journal 
is used. When the paper bears the names of more 
than one author, the journal reference appears with 
each author’s name, but the title of the paper is given 
with the senior author only. Multiple authorship is 
denoted by the symbol @ preceding the surname. 

Title: Given exactly as in the reviewing journal. 
Titles of separately bound publications (books, re- 


ports, etc.) are italicized, and are followed by the | 


publisher. 

Reference to literature: In many cases the abbrevia- 
tions of journal names have been necessarily severe 
and unconventional. The name of the journal in 
italics and the number of the volume in boldface, 
are followed by the initial page number. 
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1046 Selected Abstracts | Surface-wave resonance effect in a reactive cylin- 
y time ; ; drical structure excited by an axial line source, 
‘hittle Use of Chebychev polynomials in thin film com- | A. L. Cullen, J. Research NBS 64D, No. 1, 13 
putations, K. D. Mielenz, J. Research NBS 638A, | (1960). 


, Vv “4 

5, 811) No. 3, 297 (1909). it is shown that a purely reactive cylinder excited by a 
neighboring line source can, under suitable conditions, give 
rise to a radiation pattern closely approximating the fune- 
tion cos n @. 

In a numerical example, a cylinder of three transverse electro- 
magnetic wavelengths circumference has a surface reactance 
chosen to emphasize the term cos 6 @ in the Fourier series of 
the resultant radiation pattern. It is shown that only 1.1 
percent of the total power delivered to the line source is 
radiated in unwanted modes. 

It is also shown that the position of the line source does not 
affect this result to first order provided that k(b—a)< <1, 
where b—a is the distance of the line source from the cylindri- 
cal surface. 


Publications of the National Bureau of Standards* 


(Including Papers in Outside Journals) 


‘tuels, from Herpin’s expression for the mth power of a multilayer 
Parig matrix, very simple closed formulas are derived for the 
l, 258 matrices and optical constants of any multilayer with a 
; y periodic structure, 
rages, ‘\ecording to Epstein’s theorem, any symmetrical multilayer 

is equivalent to a fictitious monolayer. A simple expression 
» 674 for the equivalent index and thickness of this monolayer is 
auto. , deduced for the case of a periodic and symmetrical sequence 
lei of equally thick films. ; 
Ns, As compared to any other method of numerical computation, 
» 261 the suggested formulation provides a considerable saving 
arias. of time and work. In a numerical example, this saving 

9g amounts to about 80 percent. 


oeg Theory of formation of polymer crystals with | An exact earth-flattening procedure in propagation 
» 368 folded chains in dilute solution, J. I. Lauritzen, Jr., around a sphere, B. Y. -C. Koo and M. Katzin, 
dskr. snd J. D. Hoffman, J. Research NES 64A, No. 1, | J- Research NBS 64D, No. 1, 61 (1960). 


73 (1960). By a refinement of the procedure used in the usual earth- 
flattening approximation, the problem of propagation around 
a spherical earth is reduced to an exact equation of the same 
form. Thereby the earth-flattening procedure becomes appli- 
cable to arbitrarily large heights and distances. It is also 
found that existing solutions of the approximate equations can 
be reevaluated to yield the exact solutions for slightly different 
refractive index distributions. 


atist. : : ' . : 
A detailed interpretation of the kinetics of homogeneous 


. nucleation and growth of a polymer from dilute solutions 

/Sta- isgiven. The probability of forming both nuclei with folded 

chains, and the conventional bundle-like nuclei, from dilute 

atist, | solution is carefully analyzed. It is predicted that at suffi- 

ciently high dilution, critical nuclei of length I* will be formed 

’ } from single polymer molecules by sharp folding of the chain 

unty backbone. After growth, the resulting crystal is flat and 

451). platelike, the loops formed by the chain folds being on the | Some applications of statistical sampling methods 

. 167 upper and lower surfaces. Kinetic factors determine that to outgoing letter mail characteristics, N. C. Severo, 

the distance between the flat surfaces in the grown crystal . . , M. Zel VRS 

. will vary over only a narrow range of values about a charac- | A. K. Newman, S. M. oung, and Mi. Geren, 4 S 
y of teristic step height, 1*=o,/Af. Here o, is the free energy | Tech. Note 16 (PB1513875) $2.75. 

required to form a unit area of the flat surface, and possesses . 

347 a-rather large value because the work involved in forming 


This paper presents applications of statistical sampling pro- 
cedures especially devised to procure information about the 


Roy. folds is included. The quantity Af is the free energy differ- characteristics of outgoing letter mail. The results of four 
“* ence per unit volume of erystal between the erystalline and eparate studies carried out in the Washington, D.C., San 

. 7 ‘a as Tage ‘ : ” es . 42 , : oC ara > b ‘ : re v , . 7] a . 
sta-| = Sp ade — ” Ap These ly . Lt ppemcagnone 8 od hem | Francisco, and Los Angeles post offices are herein summarized. 
‘aa a . he re ry — aad sgt ay ragga sie ome ‘ly ape: | The techniques used in the various studies were developed so 
Sei ie fogarithim of the nucleation rate is approximately Pro- | that the required information would be of predetermined 


102 eo ~s (AT)? “aged “5 melting Peicted gy reliability and could be gathered without the use of a large 
ph Bs Raia To the ‘cuenta “aeons © pines og staff and without interrupting the flow of mail. The four 
the predictions given agree with the experimental results studies presented concern: (1) Letter size and color char- 
dhieised by Keller and O'Conner and ethers Gn dinate acteristics, (2) ratio of hand canceled mail to machine canceled 
etietiie of unbranched polyethylene grown frets @iate mail, (3) top and bottom clearance space of an addressed 
seletien. . . envelope, (4) proportions of long and short letters. 


23) Effect of antenna size on gain, bandwidth, and | Terrestrial propagation of very-low-frequency radio 
efficiency, K. G. Harrington, J. Research NPS 64D, | waves, J. R. Wait, J. Research NES 64D, No. 2, 153 


No. 1,1 (1960). | (1960). 

A theoretical analysis is made of the effect of antenna size | A self-contained treatment of the wave-guide mode theory 
oh parameters such as gain, bandwidth, and efficiency. | of the propagation of very-low-frequency radio waves is 
Both near-zone and far-zone directive gains are considered. | presented. The model of a flat earth with a sharply bounded 
It is found that the maximum gain obtainable from a broad- | and homogeneous ionosphere is treated for both vertical and 


band antenna is approximately equal to that of the uniformly | horizontal dipole excitation. The properties of the modes 
illuminated aperture. If higher gain is desired, the antenna | are discussed in considerable detail. 

must necessarily be a narrow-band device. In faet, the | The influence of earth curvature is also considered by refor- 
input impedance becomes frequency sensitive so rapidly that, | mulating the problem using spherical wave functions of 
for large antennas, no significant increase in gain over that complex order. The modes in such a curved guide are investi- 
of the uniformly illuminated aperture is possible. Also, if | gated and despite the initial complexity of the general solution, 
the antenna is lossy, the efficiency falls rapidly as the gain | many interesting and limiting cases may be treated in simple 
is increased over that of the uniformly illuminated aperture. | fashion to yield useful and convenient formulas for calculation. 
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Other factors considered are the influence of the earth’s 
magnetic field, antipodal effects, resonator type oscillations, 
and the influence of stratification at the lower edge of the 
ionosphere. 


Inclusion theorems for congruence subgroups, M. 
Newman and I. Reiner, Zrans. Am. Math. Soc. 91, 
No. 3, 369 (1959). 

Let G, be the group consisting of all (<¢ integral matrices of 
determinant 1. For a fixed partition ‘=r-+s of ¢ into positive 
integers r, s, and for a fixed positive integer n, define the 
subgroup 





{ ‘Aw B 


G,,.(n) = | C D® 


)eG C=0 (mod n) : 


Suppose that m is a positive integer and H a group such that 
G,,,(mn) CH CG, ,(n). 
Then 


H=G,,,(dn) where d|m. 


Further results of this type are obtained. 


New approach in the theory of satellite orbits, J. P. 
Vinti, Phys. Rev. Letters 3, No. 1, 8 (1959). 


A" gravitational potential has been found, expressed in oblate 
——— coordinates, which results in separability of the 
Hamilton-Jacobi equation for the motion of an earth satellite 
and which is much closer to the empirically accepted one than 
any heretofore used as the starting point of a calculation. 
It thus becomes possible to do orbit theory for unretarded 
satellites very accurately without the use of perturbation 
theory. 


Calculated patterns of slotted elliptic-cyiinder an- | 
tennae, J. KR. Wait and W. E. Mientka, Appl. See. 
Research B7, 449 (1959). 


The model assumed is a perfectly conducting elliptie evlinder 
of infinite length which has a narrow axial slot of finite length. 
Patterns are presented for various ratios of the major to 
minor axis. It is indicated that for reasonably large elliptic 
cylinders, the patterns are dependent mainly on the surface 
curvature in the neighborhood of the slot. 


On prediction of system behavior, J. R. Rosenblatt, 
Proc. of the N.Y. Univ—Ind. Conf. on Reliability 
Theory (New York, N.Y.) p. 89 (1988). 


It is customary in theoretical discussions of prediction of 
system performance to assume that the relation between 
variables by which performance is assessed and variables de- 
scribing the parts of the system can be specified by a function. 
The purpose of this paper is to consider some of the problems 
which arise in choosing the form of such a function and in 
specifying the relevant variables and the nature of appro- 
priate experiments, Measurements, and data. 

In the context of each of a number of familiar types of mathe- 
matical model, the possible consequences of alternative choices 
of time units, subsystem definitions, and independence 
assumptions are discussed. 


Kinetic equation for a plasma with unsteady corre- 
lations, C. M. Tchen, Phys. Rev. 114, 394 (1959). 


As a generalization of the Boltzmann equation, the kinetic 
equation for a plasma is derived in the form of a generalized 
Fokker-Planck equation, by considering unsteady correlations, 
including non-Markovian and nonlinear behavior. Both the 
binary and ternary correlations are used for many kinds of par- 
ticles with different temperatures. The coefficients of the 
kinetic equation depend on the law of interaction for a pair of 
particles and are influenced by relaxation. The effective 
potential of friction consists of two parts: the statie part 
corresponds to the Debye potential and is isotropic, the dyna- 
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mical part is axially symmetrical about the direction y 
motion, and causes a dynamical friction. The results sho, 
that the friction is proportional to velocity for slow particle 
and inversely proportional to the square of velocity for fos 
particles. This tendency of the fast particles to overeon, 
repulsion is a property connected with the “run-away” » 
electrons. <A criterion for maximum friction is derived, The 
triplet interaction, which mainly affects the shielding phe 
nomena, assures the convergence of the coefficients in case y 
distant interaction. Sinee the length scales of interactioy 
are well determined in this way, the kinetic equation can by 
expected to be valid over a longer range than does the Bolt, 
mann equation. The large scale agrees with the Debye radiys 
when the shielding term is linearized, as should be expeeted 
When time relaxation is left aside and linearization is made 
the kinetic equation degenerates to the classical Fokker.’ 
Planck equation with convergent coefficients. 
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and related functions, NBS Applied Math. Series 
50, $3.25. 
F requene vy dependence of VHF ionospheric scattering, 
Blair, NBS Tech. Note 9 (PB151368) 75 
ce Be 

Calculated behavior of a fast neutron spectrometer 
based on the total absorption principle, J. E. 
Leiss. NBS Tech. Note 10 (PB151369) $1.00. 

Penetration of gamma rays from isotropic sources 
through aluminum and concrete, M. J. Berger 
and L. V. Spencer, NBS Tech. Note 11 (PB151370) 
50 cents. 

Analysis of ionospheric vertical soundings for electron 
density profile data. Facilities for convenient 
manual reduction of ionograms, J. W. Wright 
and R. B. Norton, NBS Tech. Note 14 (PB151373) 
50 cents. 

Prediction of the cumulative distribution with time 
of ground wave and tropospheric wave trans- 


mission loss. Part 1. The prediction formula, 
P. L. Riee, A. G. Longley, and K. A. Norton. 
NBS Tech. Note 15 (PB151374) $1.50. 


Radio noise data for the International Geophysical 


Year July 1, 1957—December 31, 1958, W. Q. 
Crichlow, C, A. Samson, R. by Disney, and M. A. 
Jenkins, NBS Tech. Note 18 (PB151377) $2.50. 


Analysis of ionospheric vertic al Phas se for electron 
density profile data. II. Extrapolation — of 
observed electron density profiles above hmax F2, 

W. Wright. NBS Tech. Note 19 (PB151378) 
50 cents. 

Variations of gamma cassiopeiae, 5S. R. Pottasch, 
NBS Tech. Note 21 (PB151380) 75 cents. 

Precise time synchronization of widely separated 
clocks, A. H. Morgan, NBS Tech. Note 22 
(PB151381) $1.50. 

Design of single frequency filters, F. F. Fulton, Jr., 
NBS Tech. Note 23 (PB151382) 50 cents. 

Communication theory _ie of television band- 
width conservation, W. C. Coombs, NBS Tech. 
Note 25 (PBI5L 384) 50 cents. 

Aerodynamic phenomena in stellar atmospheres—A 


bibliography, NBS Tech. Note 30 (PB151389) 
$1.25. 

The changing character of chemical research in 
government, E. Wichers, The Chemist XXXVI, 


No. 7, 260 (1959). 

Study of the setting of plaster, K. D. Jorge ~~ and 
A.”S. Posner, J. Dental Research 38, No. 3, 491 
(1959). 











Radiation attenuation data, H. O. Wyckoff, Radia- | 
tion Hygiene Handbook, edited by H. Blatz, 
Section 8, p. 1 (McGraw-Hill Book Co., New 
York, N.Y., 1959). 


Corrosion of type 310 stainless steel by synthetic 
fuel oil ash, H. L. Logan, Corrosion 15, 443t 
(1959). 

Physical research, Part 2, G. C. Paffenbarger and 
W. Souder, J. Am. Dental Assoc. 58, 98 (1959). 
Effects of gamma radiation on collagen, J. Cassel, 
J. Am. Leather Chemists’ Assoc. LIV, No. 8, 

432 (1959). 

Standards for neutron flux measurement and neutron 
dosimetry, R. S. Caswell, E. R. Mosburg, Jr., and 
J. Chin, Dd United Nations International Conf. on 


the Peaceful Uses of Atomic Energy, Vol. 21, 
Health and Safety: Dosimetry and Standards, 
P/752 USA, p. 92 (1959). 


Absorption of radiation by a cylindrical sample of 


a strong absorber, P. H. Fang and I. A. Stegun, 
Letter to Editor, J. Chem. Phys. 31, 267 (1959). 

PILOT—A_ new multiple computer system, A. L. 
Leiner, 6 A. Notz, J. L. Smith, and A. Wein- 
berger, J J. Assoc. Computing Mach. 6, No. 3, 313 
(1959). 


Properties of rutile (titanium dioxide), F. A. Grant, 
Rev. Modern Phys. 31, 646 (1959). 

Analysis of vibrational relaxation data in shock wave 
experiments, K. E. Shuler, J. Chem. Phys. 30, 
No. 6, 1631 (1959). 

The interlaboratory ee of testing methods, 
J. Mandel and T. W. Lashof, ASTM Bull. 
239, 53 (1959). 

Branched-chain higher sugars. I. A 9-aldo-4-C- 
formyl-nonose derivative, R. Schaffer and H. S 
Isbell, J. Am. Chem. Soc. 81, 2178 (1959). 

The structure of electrolytic solutions, W. J. Hamer 
(John Wiley & Sons, Inc., New York, N.Y., 
1959). 

The use of a vacuum microbalance in studies of 
electron tube materials, G. F. Rouse, Proce. 4th 
Natl. Conf. on Tube Techniques, September 10-12, 
1958, sponsored by the Advisory Group on Elec- 
tron Tubes, p. 262 (New York University Press, 
New York, N.Y., 1959). 

Widths of cracks in concrete at the surface of rein- 
forcing steel evaluated by means of tensile bond 
specimens, D. Watstein and R. G. Mathey, J. Am. 
Concrete Inst. 31, No. 1, 47 (1959). 

Exposure standards and radiation protection regula- 
tions, L. S. Taylor, Radiation Hygiene Handbook, 
edited by H. Blatz, Section 3, p. 2 (MleGraw-Hill 
Book Co., New York, N.Y., 1959). 


Geometrical anisotropy of magnetic materials in 
wave guides and cavities, A. L. Steinert, J. Appl. 


Phys. 30, 1109 (1959). 
Radioactivity standardization in the United States, 


W. B. Mann and H. H. Seliger, 2d U.N. Intern. 
Conf. on the Peaceful Uses of Atomic Energy, 
Vol. 21, Health and Safety: Dosimetery and 


Standards, P/750 USA, p. 90 (1959). 

Magnetic interaction of Hg, V. Griffing J. L. Jackson, 
and B. J. Ransil, J. Chem. Phys. 30, No. 4, 1066 
(1959). 
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Program of the International Commission on radio- 
logical units and measurements, L. S. Taylor, 
L. H. Gray, and H. O. Wyckoff, 2d U.N. Intern, 


Conf. on the Peaceful Uses of Atomic Energy 
gy, 


Vol. 21, Health and Safety: Dosimetry and 
Standards, P/2243 WHO, p. 81 (1959). 
The nature, cause and effect of the porosity jn 


Influence of gas bubbles on 
Ogburn and D. W. 


electro-deposits, lV. 
the formation of pores, F. 


Ernst, Plating 46, 957 (1959). 

Low even configurations in the first spectrum of 
ruthenium (Rui), R. E. Trees, J. Opt. Soe. Am, 
49, 838 (1959). 

End plate modification of X-band TE), cavity 
resonators, M. C. Thompson, Jr., F. E. Freethey 
and D. M. Waters, Letter to editor, IRE Trans, 


on Microwave Theory Tech. MTT=7, 388 (July 
1959). 
Introductory remarks, A. V. Astin, Proe. 5th Tech, 
Session on Bone Char 1957, p. 1 (Bone Char 
Research Project, Inc., Charlestown, Mass., 1959), 
Analysis of liquid sugars, E. J. McDonald, Proe, 
5th Tech. Session on Bone Char 1957, p. 77 (Bone 
Char Research Project, Inc., Charlestown, Mass., 


1959). 
Development of a new test for the abrasion hardness 


of bone char, F. G. Carpenter, Proe. 5th Tech. 
Session on Bone Char 1957, p. 99 (Bone Char 
Research Project, Inc., Charlestown, Mass., 1959), 


Effects of controlled decarbonization on the per- 
formance of service synthad, F. W. Schwer, W. V. 
Loebenstein, and E. P. Barret, Proce. 5th Teeh. 
Session on Bone Char 1957, p. 145 (Bone Char 
Research Project, Inc., Charlestown, Mass., 1959). 

Some mechanisms of color and ash removal by bone 
char, A. Gee, Proc. 5th Tech. Session on Bone 
Char 1957, p. 163 (Bone Char Research Project, 
Inc., Charlestown, Mass., 1959). 

Survey of variations in use for conducting laboratory- 
scale column filteration tests, V. R. Deitz, Proe. 
5th Tech. on Bone Char 1957, p. 237 
(Bone Char Research Project, Inc., Charlestown, 
Mass., 1959). 

Comparison of column decolorization experiments 
with theory, W. V. Loebenstein, Proc. 5th Tech 
Session on Bone Char 1957, p. 253 (Bone Char 
Research Project, Inc., Charlestown, Mass., 1959). 

Sugar retention by char, F. G. Carpenter, Proce. 5th 
Tech. Session on Bone Char 1957, p. 279 (Bone 
Char Research Project, Inc., Charlestown, Mass., 
1959). 

Removal of organic anions by bone char, V. R. Ditz 
and H. M. Rootare, Proc. 5th Tech. Session on 
Bone Char 1957, p. 297 (Bone Char Research 
Project, Ine., Charlestown, Mass., 1959). 

Recent advances in cryogenic engineering, R. B. 
Jacobs, Am. Rocket Soc. ARS J. 29, 245 (1959). 


Session 


Transistorized velocimeter for measuring speed in 


the sea, C. E. Tschiegg, J. Acoust. Soc. Am. 31, 


No. 7, 1038 (1959). 

Magnetic study of the frozen products from the 
nitrogen microwave rary” B. J. Fontana, J. 
Chem. Phys. 31, No. 1, 148 (1959). 
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The thermal E. M. F. of several thermometric 
alloys, R. L. Powell and M. D. Bunch, Suppl. 
Bull. Inst. Intern. du Froid (Delft, Holland) 
Comm. 1, 129 (1958). . . 

Electrical discharge induced luminescence of solids 
at low temperatures, L. J. Schoen and R. E. 
Rebbert, J. Mol. Spect. 3, No. 4, 417 (1959). 

The nature, cause and effect of the porosity in 
electrodeposits. V. An evaluation of the sensi- 


tivity of the ferroxyl test, F. Ogburn, D. W. Ernst, | 


and W. H. Roberts, Plating 46, 1052 (1959). 

Spectrum of thin target bremsstrahlung bounded by 
a forward circular cone, J. H. Hubbell, J. Appl. 
Phys. 30, No. 7, 981 (1959). 

Negative atomic ions, H. R. Johnson and F. Rohrlich, 
J. Chem. Phys. 30, 1068 (1959). 

Factors affecting modulation techniques for VHF 





scatter systems, J. W. Koch, IRE Trans. Commun. | 


Syst. CS-7, No. 2, 77 (1959). 

The very-low-frequency emissions generated in the 
earth’s exosphere, R. M. Gallet, Proc. IRE 4%, 
No. 2, 211 (1959). 

Reaction of hydrogen atoms with solid oxygen at 
20° K., R. Klein and M. D. Scheer, J. Chem. 
Phys. 31, No. 1, 278 (1959). 

Application of the Williams-Landel-Ferry equation 
to silicate glasses, A. B. Bestul, Glastech. Ber. 
(Frankfurt, Germany) 32 K, No. VI, 59 (1959). 

Glow discharge spectra of copper and indium above 
aqueous solutions, D. E. Couch and A. Brenner, 
Tech. Notes, J. Electrochem. Soc. 106, 628 (1959). 

Free radical chemistry, J. W. Moyer and A. M. 
Bass, Chem. and Eng. News, p. 51 (Aug. 24, 1959). 

Gaseous heat conduction at low pressures and 
temperatures, R. J. Corrucini, Vacuum 7, 8, 19 
(Pergamon Press Ltd., London, England, 1959). 

Discussion of the papers of Messrs. Satterthwaite 
and Budne, W. J. 
2, 157 (1959). 

Production of embossing plates from texture patterns 
by electroforming methods, J. P. Young and V. A. 
Lamb, Plating 46, 1033 (1959). 

Method of evaluating the clinical effect of warping 
a denture: Report of a case, J. B. Woelfel and 
G. C. Paffenbarger, J. Am. Dental Assoc. 59, 25U 
(1959). 

The sound transmission loss of some 
construction, R. V. Waterhouse, R. D. Berendt, 
and R. K. Cook, Noise Control 5, No. 4, 40 (1959). 

Cryogenic insulation, R. H. Kropschot, Am. Soe. 
Heating, Refrig. and Air-Conditioning Engrs. 1, 
No. 9, 48 (1959). 

End plate modification of X-band TEQ11 cavity 
resonators, M.S. Thompson, F. E. Freethey, and 
D. M. Waters, IRE Trans. on Microw. Theory 
Tech. MT'T-7%, 388 (1959). 

Evolution of amplified waves leading to transition 
in a boundary layer with zero pressure gradient, 
P.S. Klebanoff and K. D. Tidstrom, NASA Tech. 
Note D-195, 1 (1959). 

Microwave reflectometer, G. F. Engen and R. W. 
Beatty, IRE Trans. on Microw. Theory Tech. 
MTT-7, 351 (1959). 
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Multiple ionization of sodium vapor by electron im- 
pact, V. H. Dibeler and R. M. Reese, J. Chem. 
Phys. 31, No. 1, 282 (1959). 

Penetration and diffusion of X-rays, U. Fano, L. V. 
Spencer, and M. J. Berger, Encyclopedia of Phys. 
38, No. 2, 660 (1959). 

Problems of the experimenter, W. J. Youden, Natl. 
Conv. Trans. Am. Soc. Quality Control, p. 41 
(1959). 

Vibrational intensity distributions in the nitrogen 
afterglow, U. H. Kurzweg and H. P. Broida, J. 
Mol. Spec. 3, No. 4, 388 (1959). 

Some characteristics of VLF propagation using at- 
mospheric waveforms, W. L. Taylor and L. J. 
Lange, Recent Advances in Atmos. Elec. (Proe. 
2d Conf. on Atmos. Elec.), Portsmouth, N.H. 
May 20-23, 1958, p. 609 (Pergamon Press, Ine., 
New York, N.Y., 1958). 

A new aid for the rapid determination of absorption 
corrections by Albrecht’s method, D. K. Smith, 
Acta Cryst. (Copenhagen, Denmark) 12, pt. 6 
479 (1959). 

Measurements made by matching with known stand- 
ards, W. J. Youden, W. S. Connor, and N. C. 
Severo, Technometrics 1, No. 3, 101 (1959). 

Influence of crystallographic orientation on the 
pitting of iron in distilled water, J. Kruger, J. 
Electrochem. Soc. 106, No. 8, 736 (1959). 

Evacuated powder insulation for low temperatures, 
M. M. Fulk, Progress in Cryogenics—1, 65 (Hey- 
wood and Co., Ltd., London, England, 1959). 

A direct-reading viscometer, M. R. Shafer, Instru- 
ments and Control Systems 32, No. 7, 1044 (1959). 

Atomic weights, E. Wichers, ch. 6, Treatise on Anal. 
Chem., pt. 1, Theory and practice, vol. 1, sec. B, 
Application of chemical principles, p. 161 (The 
Interscience Encyclopedia, Inc., New York, N.Y., 
1959). 

Concept and determination of pH, R. G. Bates, ch. 
10, Treatise on Anal. Chem., pt. 1, Theory and 
practice, vol. 1, sec. B. Application of chemical 
principles, p. 361 (The Interscience Encyclopedia, 
Inc., New York, N.Y., 1959). ; 

Electrode potentials, R. G. Bates, ch. 9, Treatise on 
Anal. Chem., pt. 1, Theory and practice, vol. 1, 
sec. B, Application of chemical principles, p. 319 
(The Interscience Encyclopedia, Inc., New York, 
N.Y., 1959). 

Principles and methods of sampling, W. W. Walton 
and J. 1. Hoffman, ch. 4, Treatise on Anal. Chem., 
pt. 1, Theory and practice, vol. 1, sec. A, Analyti- 
cal chemistry: Its objectives, functions and limi- 
tations, p. 67 (The Interscience Encyclopedia, Inc., 
New York, N.Y., 1959). 

System loss in radio-wave propagation, K. A. 
Norton, Letter to Editor, Proc. IRE 47, No. 9 
1661 (1959). 

The oxide films formed on copper single crystal sur- 
faces in pure water. 1. Nature of the films formed 
at room temperature, J. Kruger, J. Electrochem. 
Soc. 106, No. 10, 847 (1959). 

Strain gauge calibration device for extreme temper- 
atures, R. M. McClintock, Rev. Sei. Instr. 30, 
No. 8, 715 (1959). 
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Variation of the thermodynamic ideal temperature 
in the polystyrene-cyclohexane system, D. Me- 
Intyre, J. H. O’Mara, and B. C. Konouck, 
Am. Chem. Soc. 81, 3498 (1959). 

Rapid insertion device for coaxial attenuators, 
A. Y. Rumfelt and R. J. Como, Rev. Sci. Instr. 
30, No. 8, 687 (1959). 

Energy dissipation in standing waves in rec oe 
basins, G. H. Keulegan, J. “Fluid Mech. » pt. 
35 (1959). 

Trapped energetic radicals, J. L. Franklin and H. P. 
Broida, Ann. Rev. Phys. Chem. 10, 145 (1959). 
Electrode potentials, p. 319, and Concept and 
determination of pH, p. 361, R. G. Bates, pt. 1, 
vol. 1, Treatise on Anal. Chem., edited by I. M. 
Kolthoff and P. J. Elving (Interscience Pub. Co., 

New York, N.Y., 1959). 

Total photoelectric cross sections of copper, molyb- 
denum, silver, tantalum, and gold at 662 kev, 
W.F. Titus, Phys. Rev. 115, No. 2, 351 (1959). 

Diode in feedback loop makes stable transistor bias 
a, G. F. Montgomery, Electronic 
7, 54 (1959). 
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Design | 


Strength of synthetic single crystal sapphire and 


ruby as a function of temperature and orie —— 
J. B. Wachtman and L. H. Maxwell, J. Am. 
Ceram. Soe. 42, No. 9, 432 (1959). 


PILOT, the NBS multicomputer system, A. L. | 
Leiner, W. A. Nots, J. L. Smith, and A. Wein- | 
berger, Proc. Eastern Joint Computer Conf. | 


(Philadelphia, Pa.) p. 164 (1958). 

Investigation of creep behavior of structural joints 
under cyclic loads and temperatures, L. Mordfin, 
N. Halsey, and G. E. Greene, NASA Tech. Note 
D. 181, p. 37 (1959). 

Dimensional stability of impregnated sole leather, 
B. H. Fouquet, J. Am. Leather Chemists’ Assoc. 
LIV, No. 10, 544 (1959). 

Theoretical problems of reliability measurement 
and prediction, J. R. Rosenblatt, Proc. of the 
Third Exploratory Conf. on Missile Model Design 
for Reliability Prediction at White Sands Missile 
Range (New Mexico) p. 29 (1959). 

Accuracy and precision: evaluation and interpreta- 
tion of analytical data, W. J. Youden, ch. 3, 
Treatise on Anal. Chem., pt. 1. Theory and 
Practice, vol. 1, sec. A, Analytical chemistry: its 
objectives, functions and limitations, p. 47 (The 
Interscience Encyclopedia, Inc., New York, N.Y., 
1959). 

Dimensional changes in systems of fibrous ma- 
cromolecules: Polyethylene, L. Mandelkern, D. E. 
Roberts, A. F. Diorio, and A. S. Posner, J. Am. 
Chem. Soc. 81, 4148 (1959). 

Microwave reflectometer techniques, G. F. Engen 
and R. W. Beatty, IRE Trans. on Microw. 
Theory and Tech. MTT-—7, No. 3, 346 (1959). 

Attenuation of seattered cesium—137 gamma rays, 
F. S. Frantz, Jr., and H. O. Wyckoff, Radiology 
73, No. 2, 263 (1959). 

Observations of the ionosphere over the south 
geographic pole, R. W. Knecht, J. Geophys. 
Research 64, No. 9, 1243 (1959). 
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Dependence of mechanical strength of brittle poly. 
crystaline specimens on emg and g~* size, 
F. P. Knudsen, J. Am. Ceram. Soe. 42, No. 8 
376 (1959). 

ony of dielectric relaxation in molecul: ar crystals, 

D. Hoffman, Archives des Sciences, Soe. Phys.! 


Vist Nat. Geneve 12, 36 (1959). 

Denture reliners direct, hard, wk resin, G 
M. Brauer, E. E. White, Jr., C. L. Burns, ‘and 
J. B. Woelfel, J. Am. Dental ai 59, 270 (1959), 


Photoproton cross sections of carbon, Ne Penner and 
J. KE. Leiss, Phys. Rev. 114, No. 4, 1101 (1959), 

Cryogenic engineering of hydroge n bubble chambers 
B. W. Birmingham, D. B. Chelton, D. B. Mann, 
and H. P. Hernandez, ASTM Bull. No. 240, 34 
(TP164) 1959. 

A porcelain reference electrode conductive to sodium 
ions for use in molten salt systems, R. J. Labrie 
and V. A. Lamb, J. Electrochem. Soc. 106, No. 10, 
895 (1959). 

Memorandum on a procedure for obtaining specteal| 
radiant intensities of tungsten-filament lamps, 
400-700 my, L. E. Barbrow, J. Opt. Soe. Am. 
Letter to Editor 49, No. 11, 1122 (1959). 

Paramagnetic resonance spectra of active species, 


Blue material from hydrazoic acid, W. B. Gager 
and F. O. Rice, J. Chem. Phys. 31, No. 2, 564 
(1959). 

Ceramic wafer tubes for modular units, C. P. 


Marsden, Section on Components and Materials, 
Electronics 32, No. 39, 94 (1959). 

Magnified and squared VSWR responses for micro- | 
wave reflection coefficient measurements, R. W. 
Beatty, IRE Trans. on Microw. Theory and Teeh. , 
MTT-7, No. 3, 346 (1959). 

Argon deposition on a 4.2° K surface, 
KF. A. Mauer, and L. H. Bolz, J. Chem. 
No. 2, 546 (1959). 


S. N. Foner, 

Phys. 31, 

Low-temperature tensile properties of copper and 
four bronzes, R. M. MeClintock, D. A. Van Gundy, 
and R. H. Kropschot, ASTM Bull, No. 240, 34 , 
(TP177) 1959. 

A nomenclature for conformations of pyranoid sugars 
and derivatives, H. S. Isbell and R. S. Tipson, 
Science 130, No. 3378, 793 (1959). 

On the convergence of the Rayleigh quotient itera- 
tion for the computation of the characteristic roots 
and vectors. III (General Rayleigh quotient and 
characteristic roots with linear elementary divi- 
sors) A. M. Ostrowski, Arch. Rat. Mech. Anal. 3, 
325 (1959). 

On the convergence of the Rayleigh quotient itera- 
tion for the computation of the characteristic roots 
and vectors. LV (Generalized Rayleigh quotient 
for nonlinear elementary divisors) A. M. Ostrowski, 
IV Arch. Rat. Mech. Anal. 3, 341 (1959). 


*Publications for which a price is indicated (except for NBS 
Technical Notes) are available only from the Superintendent of 
Documents, U.S. Government Printing O flice, Washington 25, D.C. 
(foreiqn postage, one-fourth additional). Technical Notes are 
available only from the Office of Technical Services, U.S. Depart- 
ment of Commerce, Washington 25, D.C. (Order by PB number). 
Reprints from outside journals and the NBS Journal of Research 
may often be obtained directly from the authors. 


82 
O 





poly. 
Size, 
0. § 


tals, 


*hys, 


i, G, 
and 
159). 
and 
159), 





bers, ’ 


ann, 
), 34 


lium 
brie 
. 10, 


tral 
nps, 


Am. 


cles, 
ger 
564 
P 
lals, 
cro- 
W. 
ech. 


ner, 


31, 


and 
dy, 


- 


34 


rars 
on, 


sTa- 
Ots 
ind 
ivi- 


. 3, 


Ta- 
ots 
ent 


“ 


st] 


"BS 
t of 
dC. 
are 
art- 
er). 
rch 








5o- 1 6001 4 + (é 








AVCHIGAN #f y 

L AB * \e -* 
UNITE! ere MICH one ee PAYMENT OF POSTAGE, somt™t 
GOVERNMENT |} (GPO) \ 

DIVISION OF PUBLIC DOW -0is 
WASHINGTON 25, D.C. 

OFFICIAL BUSINESS 
JOURNAL OF RESEARCH of the National Bureau of Standards Vol. 64B, No. 1 January-March 1960 | 


B. Mathematics and Mathematical Physics 





Contents 
Page 

Theory of diffraction in microwave interferometry. D. M. Kerns and 

E. S. Dayhoff. ] 

Some solutions for electromagnetic problems involving spheroidal, spheri- 

cal, and cylindrical bodies. James R. Wait. 3s @ 

Kantorovich’s inequality. Morris Newman. 33°47 

A symmetric continuous poker model. A. J. Goldman and J. J. Stone. 35 

Moebius function on the lattice of dense subgraphs. R. E. Nettleton and 

M. S. Green. 4] 

The minimum of a certain linear form. Karl Goldberg. 49 

Space of k-commutative matrices. Marvin Marcus and N. A. Khan. 51 

Selected bibliography of statistical literature, 1930 to 1957: I. Correla- | 

tion and regression theory. Lola S. Deming. 5S | 

Selected bibliography of statistical literature, 1930 to 1957: II. Time - 

series. Lola S. Deming. 69 

Publications of the National Bureau of Standards. 77 
ES SET 





For sale by the Superintendent of Documents, U.S. Government Printing Office, Washington 25, D.C. Price 75 
cents (single copy). Subscription price: $2.25 a year, $0.50 additional for foreign mailing 


UNITED STATES GOVERNMENT PRINTING OFFICE, WASHINGTON : 1960 





